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SEEFAR:  AN  IMPROVED  MODEL  FOR  PRODUCING  LINE-OF-OF-SIGHT-MAPS 


1.  INTRODUCTION 

A  1 ina-of-sight  map,  or  intervisibility  plot,  is  a  graphical  representa¬ 
tion  of  those  portions  of  a  geographical  area  which  are  visible  to  a  given 
observer.  (The  term  "visible"  throughout  this  report  refers  to  the  existence 
of  line  of  sight  from  an  observation  point  to  a  target  point  which  is  uninter¬ 
rupted  by  intervening  terrain  or  obstacles.)  These  maps  may  be  as  shown  in 
Figure  1,  where  the  hatched  portions  represent  those  areas  in  which  a  target 
of  a  given  height  would  be  obscured.  Targets  in  the  blank  portions  would  be  . 
visible.  Or  the  maps  may  be  shown  as  in  Figure  2,  where  multiple  target 
heigfits  are  considered.  Each  symbol  represents  a  particular  target  height. 

A  symbol  plotted  in  a  certain  area  indicates  that  targets  of  the  associated 
height  (or  higher)  located  in  this  area  are  visible  to  the  observer.  For  ex¬ 
ample,  in  Figure  2  a  target  10  meters  high  (or  higher)  would  be  visible  to  the 
indicated  observer  if  located  at  the  point  with  UTM  (Universal  Transverse 
Mercalor  Projection)  coordinates  (554500,  5612600). 

A  number  of  computer  programs  are  available  for  providing  these  maps; 
most  of  these  programs  are  based  on  the  algorithm  described  as  "an  intuitive 
approach"  in  this  report,  LOSMAP  is  one  such  program  available  at  USAMSAA. 

Use  of  an  alternate  algorithm,  however,  has  resulted  in  a  considerable  re¬ 
duction  in  both  the  memory  requirements  and  compute  time.  An  implementation 
of  this  algorithm,  called  SEEFAR,  is  included  in  Appendix  A.  The  purpose  of 
this  report  is  to  describe  the  basic  SEEFAR  algorithm  and  to  compare  it  to 
the  LOSMAP  approach. 

Line-of-sight  maps  have  proved  particularly  useful  in  the  development  of 
combat  scenarios  for  selecting  reasonable  observer  positions  and  tactically 
sound  attack  routes.  These  maps  can  be  useful  to  the  wargamer,  to  the  tac¬ 
tician  in  the  field,  and  to  the  test  plan  developer.  Additionally,  they  may 
be  helpful  in  developing  desirable  weapon  system  characteristics.  For  example, 
it  would  be  fruitless  to  develop  a  weapon  system  requiring  line  of  sight  with 
a  range  of  10  kilometers  to  be  used  in  an  area  where  observers  typically  can¬ 
not  see  a  target  at  further  than  5  kilometers. 

The  information  required  to  produce  these  maps  is  of  three  types.  First, 
the  positional  relationship  between  the  map  and  observer  must  be  specified. 

An  area  of  interest  is  denoted  by  its  extreme  rectangular  coordinates.  The 
observer's  position  must  be  within  this  area.  Secondly,  heights  of  the  ob¬ 
server  and  the  target  are  required.  The  target  height  is  the  height  of  that 
point  on  the  target  that  must  be  seen  before  the  target  is  considered  visible; 
the  observer  height  designates  the  eye  or  sensor  position  above  the  ground. 
Finally,  the  actual  terrain  elevations  along  with  any  ancillary  terrain  data 
must  be  supplied. 

The  Defense  Mapping  Agency  (DMA)  digitized  terrain  data  fulfill  the  third 
information  requirement.  These  data  were  created  by  extracting  elevations 
from  concour  lines  on  1:50000  scale  contour  maps.  Then,  through  planar  inter¬ 
polation,  elevations  were  obtained  at  12.5-meter  intervals.  This  information 
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Figure  1.  Line -of -Sight  Map. 
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is  supplied  on  magnetic  tapes  as  south- to-north  strings  of  data  sweeping  from 
west  to  east.  The  term  scanline,  or  scan,  throughout  this  paper  will  denote 
a  south- to-north  string  of  information.  Thus,  the  DMA  data  is  a  set  of  eleva¬ 
tion  scanlines.  Additional ly>  an  indication  of  the  presence  of  a  forest, 
orchard  or  urban  area  has  been  included  with  each  elevation. 
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2.  AN  INTUITIVE  LINE-OF-SIGHT  ALGORITHM  (LOSMAP) 


2.1  Algorithm  Description. 

Now,  given  the- information  described  in  the  previous  section,  suppose 
we  wish  to  produce  a  map  like  that  in  Figure  1.  Basically,  at  evenly  spaced 
intervals  thoughout  the  area  of  interest  the  question  must  be  asked,  "Can  the 
observer  see  the  target  if  it  is  at  this  spot?"  A  reasonable  approach  to 
answering  this  question  is  to  draw  the  terrain  profile  between  the  observer  and 
the  target,  then  draw  the  line  from  the  observer  eyeball  to  the  target  point 
and  determine  whether  the  terrain  profile  interrupts  that  observer-target  line. 
If  there  is  no  interruption,  line  of  sight  exists;  otherwise,  this  spot  is  not 
visible. 


In  Figure  3  the  observer  and  target  are  depicted  in  the  array  of 
digitized  elevations.  As  the  figure  indicates,  at  each  point  where  the  line 
between  observer  and  target  crosses  a  grid  line,  linear  interpolation  is  used 
to  find  another  height  on  the  desired  profile.  As  the  profile  is  thus  being 
constructed,  a  comparison  is  made  between  the  elevation  of  each  new  point  on 
the  profile  and  the  corresponding  height  of  the  line  drawn  from  the  observer 
eyeball  to  the  target  point.  If  the  profile  height  is  ever  the  larger  of  the 
two  heights,  line  of  sight  does  not  exist;  otherwise,  line  of  sight  does  exit 


Of  course  this  is  a  simplification  of  the  problem.  Because  the  pre¬ 
sence  of  vegetation  and  urban  areas  can  have  a  major  impact  on  visibility  for 
an  observer,  it  is  generally  desirable  to  use  the  DMA  vegetation/urban  in¬ 
dicators.  For  example,  while  constructing  the  profile  an  appropriate  height 
increment  may  be  added  to  the  profile  if  it  runs  through  trees  or  towns.  Fur¬ 
ther,  the  effect  of  earth  curvature  must  be  considered.  If  multiple  target 
heights  are  considered,  some  accounting  method  must  be  used  to  keep  track  of 
the  minimum  visible  target  height.  This  brief  description  gives  enough  infor¬ 
mation  to  point  out  the  major  problems  with  the  algorithm:  computer  space 
and  time  requirements. 


2 . 2  Problems  Encountered:  Storage  and  Compute  Time 

As  long  as  the  area  of  interest  is  small,  there  is  no  problem.  The 
use  of  high  speed  computers  enables  one  to  perform  the  required  calculations 
with  ease. 

Notice,  however,  that  this  algorithm  requires  easy  access  to  elevations 
in  the  entire  area  of  interest.  For  each  target  point,  an  observer-to- target 
profile  is  constructed,  requiring  the  elevations  of  many  intervening  points. 
Consider  making  a  6  x  6  kilometer  map.  If  a  12.5-meter  grid  is  used,  an  array 
of  elevations  dimensioned  481  x  481  is  needed.  Storing  one  elevation  per  word 
would  require  231,361  words  of  memory  for  the  elevation  array  alone!  This  pro¬ 
blem  can  be  somewhat  alleviated  by  using  a  thinned  grid  of  25,  50  or  even  100 
meters.  While  this  thinning  may  result  in  more  inaccuracy,  the  outcome  is  often 
deemed  adequate.  Further,  these  elevations  may  be  packed  into  the  elevation 
array  several  at  a  time  thereby  decreasing  the  required  storage  even  more.  But, 
inevitably,  as  the  size  of  the  area  that  can  be  handled  is  increased,  a  problem 
is  encountered  that  requires  an  area  slightly  larger.  Some  have  solved  the 


storage  problem  by  swapping  data  in  and  out  of  memory  as  needed.  This,  however, 
is  rather  cumbersome  and  will  tremendously  increase  input-output  time  require¬ 
ments.  As  one  becomes  concerned  with  larger  areas  of  interest,  as  in  air-to- 
ground  and  air-to-air  intervisibility  studies,  this  storage  problem  must  be 
faced  some  other  way. 

Another  problem  to  consider  is  that  of  compute  time.  If  one  needs  a 
simple  1 ine-of-sight  map  like  that  in  Figure  1,  it  is  hard  to  predict  exactly 
how  many  calculations  will  be  performed.  For  each  point  on  the  map  a  profile 
must  be  started,  so  for  an  n  x  n  map,  at  least  n  calculations  will  be  performed. 
If  the  profile  interrupt  is  encountered  early,  few  extra  calculations  will  be 
needed.  But,  when  line  of  sight  to  the  target  point  exists,  the  entire  profile 
must  be  constructed  and  each  height  comparison  must  be  performed.  With  maps 
like  that  in  Figure  2,  considering  multiple  target  heights,  it  becomes  more 
likely  that  the  entire  profile  must  be  constructed  for  e\^ery  target  point.  The 
required  number  of  calculations  would  then  be  of  order  n  (See  Appendix  Cl). 

For  large  areas,  this  could  easily  be  a  prohibitively  expensive  algorithm. 
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3.  DYNAMIC  PROGRAMMING  ALGORITHM  (SEEFAR) 


3.1  Algorithm  Description. 

A.  One  Ray.  Now  consider  an  alternate  approach  to  the  development 
of  line-of- sight  mapL  How  should  line-of-sight  calculations  be  performed  if 
all  the  target  positions  lay  along  one  line?  For  each  target  position  along 
thf!  ray  the  question  must  be  answered,  "How  high  must  the  target  be  raised  to 
be  seen?"  But  it  hardly  seems  necessary  to  draw  a  new  profile  for  each  target 
posHion.  Instead,  one  might  choose  to  work  from  the  observer  out,  extending 
the  profile  to  each  new  target  position,  not  keeping  track  of  all  the  previ¬ 
ously  encountered  elevations  but  simply  updating  the  effect  of  the  maximum 
interrupt  encountered  thus  far  as  new  interrupts  are  encountered.  This  maxi¬ 
mum  interrupt  effect  projected  to  the  current  target  position  is  essentially 
the  minimum  target  height  required  to  overcome  the  effect  of  terrain  between 
the  observer  and  the  target.  This  height  can  be  thought  of  as  a  running 
horizon,  running  in  the  sense  that  it  must  be  updated  as  target  positions  move 
further  from  the  observer. 

For  each  target  position,  then,  the  following  procedure  could  be 

followed: 


(1)  Check  to  see  if  the  current  target  position  is  above  or  below 
the  current  horizon  (i.e.,  in  or  out  of  view) 

(2)  Plot  the  results 

(3)  Update  the  running  horizon  to  reflect  its  movement  outward,  and 
any  greater  interrupts  encountered.  (This  new  horizon  value  will  be  used  in 
the  line  of  sight  calcuation  for  the  next  target  position.) 

Consider  now  this  running  horizon,  or  maximum  interrupt  effect.  Fig¬ 
ure  4  demonstrates  maximum  interrupts,  how  they  can  change  as  the  profile  is 
extended,  and  how  their  corresponding  projected  effect  must  continually  be  up¬ 
dated.  The  maximum  interrupt  effect  must  be  updated  both  as  new  interrupts 
are  encountered  and  as  the  target  position  is  changed.  A  comparison  of  Figures 
4A  and  4B  demonstrates  the  need  to  update  the  running  horizon  as  greater  inter¬ 
rupts  are  encountered.  Considering  Figures  4B  and  4C,  it  can  be  seen  that  the 
effect  of  the  maximum  interrupt  will  vary  for  different  target  positions  even 
though  the  maximum  itself  may  remain  the  same.  One  should  further  note  that, 
as  shown  in  Figure  5,  the  maximum  interrupt  is  not  necessarily  the  terrain 
profile  with  the  highest  elevation. 

This  procedure  has  several  nice  features.  First,  because  it  is  not 
necessary  to  draw  a  new  profile  for  each  target  position,  fewer  steps  are  re¬ 
quired.  Further,  it  is  not  necessary  to  store  all  the  elevations  between  the 
observer  and  the  target  to  determine  whether  line  of  sight  exists.  The  ele¬ 
vation  of  the  current  target  position  and  the  horizon  effect  are  the  only 
■values  needed.  This  method  can  be  extended  to  produce  a  full  map,  rather  than 
a  single  ray,  greatly  reducing  the  storage  and  compute  time  problems. 
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Figure  4.  Updating  the  Running  Horizon  (The  Effect  of  the 
Maximum  Interrupt  Projected  to  the  Target  Position) 
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Figure  5.  Maximum  Interrupt  Not  Always  Highest  Elevation. 


(1)  Radial  Map.  Now  how  can  this  one-dimensional  algorithm 
be  expanded  to  two  dimensions?  One  approach  would  be  to  produce  a  radial 
line-of-sight  map,  extending  many  profiles  out  from  the  observer.  Certainly 
if  line-of-sight  Information  can  be  produced  for  one  ray  as  described  in  the 
previous  section,  it  would  be  a  simple  matter  to  project  rays  from  the  observe 
at  small  equal ly-space  angle  increments,  combine  the  results  and  produce  a 
two-dimensional  map.  With  this  method,  however,  the  line-of-sight  information 
would  be  very  closely  spaced  near  the  observer  and  more  widely  spaced  farther 
from  the  observer.  As  the  profiles  got  further  apart  the  likelihood  of  miss¬ 
ing  an  important  terrain  feature  would  increase.  Assuming,  then,  that  evenly 
spaced  data  is  more  desirable,  another  method  must  be  devised  to  produce  the 
two-dimensional  maps. 

(?.)  Rectangular  Map.  As  in  the  ray  method  the  approach  will 
be  to  work  from  the  observer  out  (i.e.,  from  the  observer  eastward,  then  from 
the  observer  westward),  keeping  track  of  a  running  horizon.  But  this  time  in¬ 
stead  of  having  a  running  horizon  point,  there  will  be  a  south-to-north  string 
of  running  horizon  points  (a  horizon  scanline)  sweeping  outward  from  the  ob¬ 
server  as  the  calculations  are  performed. 

Remember,  the  map  is  actually  a  lattice  of  target  spots.  The  goal 
is  to  determine  for  each  target  spot  the  value  that  denotes  how  high  a  target 
must  be  to  be  seen.  Let  HOWHI(i)  be  the  height  the  ith  target  must  be  raised 
to  be  seen.  If  the  target  height  is  less  than  HOWHI(i)  it  cannot  be  seen; 
otherwise  it  can.  The  following  sections  describe  a  method  for  obtaining 
describe  a  method  for  obtaining  these  H0WHI(i)'s.  In  order  to  reach  this 
goal  one  must  continually  keep  up  with  HORIZON(i),  the  effect  of  the  maximum 
interrupt  encountered  between  the  observer  and  the  current  target  position. 
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Scans  Closest  to  Observer 

Because  the  plan  is  to  work  outward  from  the  observer,  the 
first  step  is  to  perform  calculations  for  the  two  target  scanlines 
closest  to  the  observer.  In  checking  for  line  of  sight  to  the  points 
on  these  scans,  the  first  iterations  of  the  east-running  horizon  and 
west- running  horizon  must  also  be  formed.  Figure  6  will  help  illustrate 
this  procedure. 

Consider  the  point  IR  (as  in  Figure  6A).  First  determine  how 
high  a  target  must  be  to  be  seen  at  position  IR.  There  are  obviously  no 
terrain  profile  elevations  between  the  observer  and  this  point,  so  one 
might  automatically  say  HOWHI(IR)  =  0.  But  if  forests  and  urban  areas 
are  to  be  considered  a  check  must  be  made  to  determine  whether  such  an 
obstacle  is  located  at  IR.  If  so,  HOWHI(IR)  =  0  +  HVEG(IR)  where  HVEG 
is  the  height  of  the  obstacle.  If  the  target  height  is  less  than  HOWHI(IR), 
a  target  at  IR  cannot  be  seen  and  the  appropriate  symbol  must  be  plotted. 

But  more  must  be  done  at  this  position.  The  running  horizon  value  at  IR 
must  be  determined.  The  highest  elevation  (vegetation  and  urban  included) 
encountered  between  the  observer  and  IR  is  simply  the  elevation  at  IR 
(there  are  no  elevations  given  elsewhere)  plus  the  obstacle  height  at  IR. 
Thus,  HUrTZON  (IR)  =  E(1R)  +  HVEG(IR)  where  E(1R)  is  the  terrain  elevation 
at  IR.  A  similar  procedure  for  position  IL  results  in  HOWHI(IL)  = 

0  +  HVEG(IL)  and  HORIZON(IL)  =  E  (IL)  +  HVEG(IL). 

Now  consider  the  point  2R  (as  in  Figure  6B).  Determine  the 
maximum  elevation  encountered  between  the  observer  and  position  2R.  This 
is  the  horizon  value  at  (X,Y),  i.e.,  Z,  where  Z  is  the  linear  interpolation 
between  the  horizon  at  IL  and  the  horizon  at  IR  (See  Appendix  B1  for  a  dis¬ 
cussion  of  the  calculation  of  X,Y,  &  Z).  Now  project  this  Z  value  to  the 
position  2R,  obtaining  Z'(2R)  (as  described  in  Appendix  B2).  In  order  to 
be  seen,  the  target  point  must  be  higher  than  Z'(2R),  so  H0WHI(2R)  = 

Z' (2R)-E(2R) .  In  addition  the  vegetation  must  be  considered.  To  guarantee 
the  target  is  above  the  ground  H0WHI(2R)  must  be  the  maximum  of  either 
HVEG(2R)  or  Z' (2R)-E(2R) .  Once  HOWHI  is  determined  the  proper  symbol  is 
plotted  for  this  point.  Now  the  horizon  at  2R  must  be  computed.  The 
horizon  will  simply  be  the  maximum  of  (the  projection  of  the  old  horizon 
to  2R)  and  (the  current  elevation,  vegetation/urban  included). 

That  is,  H0RIZ0N(2R)  =  MAX  (Z'(2R),  E(2R)  +  HVEG(2R)  ) 

Similarly,  H0WHI(2L)  =  MAX  (Z'(2L)  -  E(2L),  HVEG(2L)  ) 

and  H0RIZ0N(2L)  =  MAX  (Z'(2L),  E(2L)  +  HVEG(2L). 

In  general,  for  the  ith  point  above  the  observer  (See  Figure  6C) 

HOWHI(i)  =  MAX  (Z'(i)  -  E  (i),  HVEG(i)) 

and  HORIZON(i)  =  MAX  (Z'(i),  E(i)  +  HVEG(i)). 

Points  south  of  the  observer  should  be  handled  in  a  manner  similar  to 
those  north  of  the  observer,  resulting  in  identical  equations  for  HOWHI (i) 
and  HORIZON(i). 
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Figure  6.  Calculation  of  HOWHl  and  HORIZON  for  Scons  Closest  to 
the  Observer. 


Scans  East  of  Observer 


Consider  scanline  R«  in  Figure  7A,  consisting  of  points 
(...OR^,  IR2.  2R2,  SRg,  ...(n-liRg,  nRg).  In  order  to  obtain  HOWHI's 
for  this  line,  only  the  horizon  values  for  scanline  R  and  the  terrain 
elevations  for  Scanline  Rg  are  needed.  Working  upward  from  the  observer, 
the  first  target  point  to'^consider  is  IRg.  Check  the  terrain  elevation 
at  IRg.  This  must  be  compared  with  the  effect  of  the  maximum  interrupt 
encountered  earlier  (the  horizon  effect).  The  intersection  of  the 
observer- target  line  with  the  line  of  known  horizon  values  closest  to 
the  target  is  at  (X,Y).  A  linear  interpolation  between  HORIZON(IR) 
and  HORIZON(OR)  will  provide  Z,  the  horizon  height  at  (X,Y).  Project  Z 
to  IRg,  obtaining  V  (IRg).  Then 

HOl.'HIdRg)  =  MAX  (Z'dRg)  -  EdRg).  HVEGdRg)) 
and  HORIZON  (IRg)  =  MAX  (Z'(1R2).  EClRg)  +  HVEGdRg)). 

Similarly,  using  Figure  7B,  the  appropriate  values  for  point 
2R2  can  be  determined. 

Figure  7C  depicts  a  variation  in  this  method.  The  observer- 
target  line  intersection  with  the  line  of  known  horizon  values  closest 
to  the  target  is  not  along  line  R,  but  rather  at  a  point  (X,Y)  between 
target  points  2R  and  ZRg.  In  this  case  the  HOWHI  and  HORIZON  equations 
are  as  follows:  ^ 

HOWHiORg)  =  MAX(Z'(3R2)  -  EORg),  HVEGORg)). 

HORIZONCSRg)  =  MAXd'SRg)  +  EORg)  +  HVE6(3Rg)). 

Thus  the  intersection  of  the  observer-target  line  with  the  closest 
vertical  horizon  scan  line  is  not  necessarily  the  intersection  of  in¬ 
terest.  The  explanation  in  Appendix  B1  for  deriving  X,  Y  and  Z  takes 
both  of  the  cases  into  account. 
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Scans  West  of  Observer 


In  general,  for  po'Jnts  east  of  the  observer 

HOWHI(i)  =  MAX  (Z'(i)  -  E(i).  HVEG(i)  ) 

HORIZON(i)  =  MAX  (Z'(i).  E{i)  +  HVcG(i)  ) 

By  simply  reversing  directions  and  starting  with  the  horizon  values 
corresponding  to  scanline  L  (calculated  while  working  with  scans  closest 
to  the  observer)  the  same  equations  are  found  to  apply  to  points  west 
of  the  observer.  Variations  in  the  computation  of  values  (X,Y,Z)  in 
determining  horizon  effects  are  all  accounted  for  in  Appendix  Bl. 
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Algorithm  Outine 

In  short,  the  map  can  be  thought  of  as  a  lattice  of  equally  spaced 
target  positions.  These  are  the  points  for  which  terrain  elevations  are 
available  and  can  be  thought  of  as  south-to-north  strings  of  target 
positions,  or  target  scanlines.  The  following  provides  a  general  outline 
of  the  flow  of  the  algorithm: 

FIND  AND  PLOT  LINE  OF  SIGHT  INFORMATION  FOR  SCANS  CLOSEST  TO  OBSERVER 

WHILE  CREATING  ORIGINAL  HORIZONS  EAST  AND  WEST  OF  OBSERVER. 

FOR  SCANS  FROM  OBSERVER  TO  SIDE  BOUNDARY  (FIRST  EASTWARD  THEN 

WESTWARD  DO: 

READ  ELEVATIONS  FOR  SCAN  OF  INTEREST 

FOR  POINTS  FROM  OBSERVER  NORTH  THEN  FROM  OBSERVER  SOUTH  DO: 

FIND  INTERSECTION  OF  OBSERVER-TARGET  LINE  WITH  HORIZON 

LINE  CLOSEST  TO  TARGET 

INTERPOLATE  BETWEEN  KNOWN  HORIZON  VALUES  TO  OBTAIN  HORIZON  VALUE 

PROJECT  HORIZON  VALUE  TO  CURRENT  TARGET  POSITION  (TAKING  INTO 

ACCOUNT  THE  EFFECT  OF  EARTH  CURVATURE) 

COMPARE  CURRENT  ELEVATION  WITH  PROJECTED  HORIZON  TO  DETERMINE 

WHETHER  LINE  OF  SIGHT  EXISTS 

UPDATE  HORIZON  VALUE 

PLOT  RESULTS 


STOP 
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3 . 2  Solution  to  Storage  and  Compute  Time  Problem. 


A.  Storage.  The  savings  in  storage  requirements  for  the  SEE- 
FAR  algorithm  are  considerable.  The  profile  algorithm  required  all  the 
elevation  data  to  be  easily  accessible,  since  many  scans  of  data  were 
needed  to  produce  each  profile;  the  newer  algorithm  needs  only  the  cur¬ 
rent  scanline  of  elevation  data  and  the  most  recently  updated  horizon 
information  to  perform  the  necessary  computations.  The  amount  of  random 
access  memory  required,  then,  is  a  function  of  the  north- to-south  di¬ 
mension  of  the  map  rather  than  the  area  of  the  map.  A  map  that  is 

100  X  50  kilometers  in  size  would  take  no  more  main  memory  than  a  map 
1  X  50  kilometers,  given  equal  grid  size.  The  original  algorithm  would 
have  taken  100  times  more  storage.  Maps  that  would  have  been  virtually 
impossible  with  the  old  algorithm  are  now  quite  possible.  But  if  only 
storage  requirements  are  reduced,  the  reader  may  wonder  whether  the 
machine  time  requirements  would  be  prohibitive  for  production  of  these 
larger  maps.  The  savings  in  computing  time  resulting  from  this  approach 
must  be  considered. 

B.  Computing  Time.  As  mentioned  earlier,  the  original  al¬ 
gorithm  required  on  the  order  of  n"^  operations  for  an  n  x  n  sized  map, 

because  each  profile  consists  of  many  points,  and  each  of  these  points 

might  have  to  be  considered  in  a  single  line-of-sight  calculation.  The 
newer  algorithm,  however,  simply  compares  the  current  target  elevation 
with  the  horizon  effect  at  the  target  point.  Thus,  a  small  fixed  number 
of  calculations  are  required  for  each  target  point  on  the  map.  This 

algorithm  is,  then,  of  the  order  n^  for  an  n  x  n  map.  For  maps  of  larger 

sizes  the  computing  time  savings  will  be  considerable.  It  can  in  many 
instances  be  the  difference  between  a  map  reasonably  produced  and  one  so 
expensive  in  computation  cost  that  its  value  tOpthe  tactician  or  anlayst 
cannot  outweigh  the  expense.  Note  that  since  n*^  target  points  are  in¬ 
volved  in  an  n  X  n  map,  order  is  the  minimum  possible  complexity.  A 
few  examples  of  these  differences  are  discussed  in  the  following  section. 


4.  COMPARISON  OF  MAPS  RESULTING  FROM  EACH  ALGORITHM 


4.1  Discussion.  The  LOSMAP  and  SEErAR  algorithms  discussed  are 
quite  different  in  nature.  While  each  uses  a  method  of  linear  inter¬ 
polation  to  get  results,  the  intuitive  method  interpolates  between  known 
elevations  to  get  elevations  along  a  profile,  then  checks  to  see  if  those 
elevations  interrupt  line  of  sight.  The  dynamic  programming  approach, 

on  the  other  hand,  keeps  up  with  the  effect  of  maximum  interrupts  en¬ 
countered  via  the  horizon  array.  A  linear  interpolation  is  made  between 
horizon  values  to  find  a  projection  of  the  horizon  to  a  specific  target 
point.  It  is  not  surprising  to  discover  the  results  of  these  methods 
are  not  always  the  same.  But,  while  the  "intuitive"  LOSMAP  method  de¬ 
scribed  uses  interpolation  and  is  therefore  not  an  exact  solution  to  the 
line  of  sight  problem,  it  would  be  disappointing  to  find  large  discrep¬ 
ancies  in  the  results  of  the  two  approaches.  One  would  certainly  hope 
the  results  of  a  new  approach  would  conform  somewhat  to  those  of  an 
intuitive  method. 

As  it  turns  out,  the  results  are  quite  similar.  In  fact,  very 
close  inspection  is  generally  required  to  distinguish  differences  in 
1 ine-of-sight  maps  with  one  target  height.  It  is  easier  to  detect  dif¬ 
ferences  when  multiple  target  heights  are  used.  The  SEEFAR  approach 
might  at  one  time  indicate  a  target  visible  that  the  LOSMAP  approach  in¬ 
dicated  nonvisible.  Another  time  the  reverse  situation  might  occur.  In 
general,  however,  the  boundaries  of  in-view  and  out-of-view  areas  do  not 
differ  widely  between  algorithms.  It  should  be  added  that,  while  the 
dynamic  programming  approach  may  not  be  the  first  to  come  to  mind  when 
tackling  the  1 ine-of-sight  problem,  it  is  a  reasonable  approach,  one 
certainly  not  counter  to  intuition.  And  the  savings  in  time  and  storage 
cannot  be  ignored. 

4.2  Sample  Maps.  The  following  pages  will  provide  the  reader  an 
opportunity  to  compare  the  results  of  the  two  algorithms. 
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MAP  OF  POINTS  VISIBIE  TO  A  SINGl£  C6SS!V£R 
oaSERVFR  COORDINATES  (SOSOOO.  5613000) 
OBSERV;  (  HEIGHT  »  2000.00 
TARC-.  HEIGHT  >  1.50 

Figure  lOA.  LOSMAP  Algorithm 

Compotu  _Tim«t  576  S«c 


5.  SUMMARY 


5.1  Advantages  of  the  SEEFAR  Algorithm.  The  SEEFAR  algorithm, 
then,  provides  a  method  for  generating  line-of-sight  maps  (and  thereby 
obtaining  intervisibility  data)  for  a  much  larger  area  or  for  a  smaller 
grid  than  could  reasonably  be  provided  using  the  LOSMAP  algorithm.  pThis 
dynamic  programming  approach  results  in  an  order  of  complexity  of  r\^  for 
an  n  X  n  map,  which  is  the  minimum  that  can  be  expected.  Further,  it 
reduces  the  storage  requirements  such  that  only  one  scanline  of  elevation 
data  need  be  in  memory  at  any  one  time.  Thus,  it  is  an  order  of  magni¬ 
tude  better  than  the  LOSMAP  algorithm  in  both  time  and  storage  requirements 

5.2  The  SEEFAR  Program.  Appendix  A1  contains  a  copy  of  the  SEEFAR 
program,  an  implementation  of  the  algorithm  discussed  in  this  report. 

It  should  be  noted  that  the  algorithm  discussion  is  concerned  with  the 
basic  approach  to  the  line-of-sight  problem.  The  SEEFAR  program  has 
several  embellishments.  First,  if  an  observer  is  in  a  forest/urban  area, 
and  the  proper  option  is  selected,  SEEFAR  will  attempt  to  move  the  ob¬ 
server  out  of  the  obstruction.  Second,  the  ratio  of  visible-to-total 
targets  in  a  selected  area  of  interest  is  given.  Third,  terrain  data  on 
as  many  as  12  DMA  tapes  can  be  merged  and  put  into  a  direct  access  file. 
Finally,  the  HOWHI  array,  which  indicates  for  every  target  position  how 
high  the  target  must  be  raised  to  be  seen  by  the  selected  observer,  can 
be  saved  on  a  permanent  file.  This  allows  one,  through  short  fast¬ 
running  programs,  to  vary  target  heights  or  merge  line  of  sight  data  for 
several  observers  obtaining  maps  which  indicate  areas  visible  to  m-out- 
of-n  observers  (1-  m  -  n).  These  m-out-of-n  maps,  or  composite  maps, 
are  done  by  a  separate  program  now,  but  could  easily  be  provided  via 

the  SEEFAR  program. 

5.3  Problems  Remaining.  There  are  line-of-sight  data  requirements 
for  which  SEEFAR  does  not  seem  a  viable  alternative.  For  example,  the 
LOSPATH  program,  in  use  at  AMSAA,  checks  along  paths  for  line  of  sight. 

To  use  SEEFAR  for  this  problem  one  must  generate  a  rectangular  map  con¬ 
taining  all  the  path  points  and  then  check  for  line  of  sight  at  appro¬ 
priate  points  within  the  rectangle.  Because  SEEFAR  builds  on  information 
generated  for  all  the  points  between  the  observer  and  the  current  target 
position,  it  does  not  seem  optimal  when  target  points  are  in  scattered 
positions  rather  than  at  evenly  spaced  intervals  within  a  rectangle  of 
interest. 

The  question  of  how  much  discrepancy  to  expect  between  the 
LOSMAP  and  SEEFAR  algorithms  has  not  been  determined  analytically.  This 
may  be  of  interest;  but  since  each  approach  provides  an  interpolated 
approximation  to  the  line-of-sight  problem,  the  value  of  such  an  analysis 
is  questionable. 
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APPENDIX  A1 
SEEFAR  PROGRAM 


pr.-ICF./.M  rEETARCMPUT/  OUTPUT/  TAPET- IMPUT/  TAPuf-OUTPUT/ 

»  T/.Pf-’/  TAPECi  TaPE''*  ''APriO/  TAPill/  TAPcl.?/ 

■"APClAi.  TaPEIu/  TAPE17,  TAPFie/  TAP^19/ 

*  T*.P'i2»  tap:3i 

*  TAPCIS) 

I 

r.AP.LAP.A  CF.OO.ME 

PF.UPMCr:  Lir.n  up  CICII*  M^.P  UCTMG  nY|;,,‘^:(;  pp.UGP.Af^MZMG  ALGOF.ITHM. 
LO;'  paTa  ('ii;:  iipuu:  .ffa?)  cam  5E  r’nRr  ■'■k  r::K  '^'t.  L..Tr''  '/liPK. 
ri.LC'.'LMEi  rf.ACTifji';  OP  (vistBiL  Ar.r.M)  f  i^'i  „  a-./.a)  :r 


r,: 

C'^AnlL 

-  OF 

ti:“::p,cst. 

TAPi^lS 

ilOLo: 

PLOT 

PATh 

TAPE  2 

•lOLPS 

McrGCD  “lfp.mZm  data 

T/.  P  L  3 

MOLPS 

'ig\m: 

DmTA  FGP.  F'lTU 

FZ  MAPS 

TaPCT 

MCLPS 

C  A  /« 0 

IIIPUT 

*APZ6 

MOLDS 

PP.IMT 

EP  OUTPU’ 

ALL  ini 

1C  nr  hF 

C  DMA 

•.r.r.AZM  u*.'‘A 

tap::  (rr.L-HEr.GE ) 

COIlfMM  GF.ir'f, 

cni:i'/iM  /:ci/  xge. :/  vqb:,  riPc:#  obtraII 

COMKHM  /AF  ’A/  XUP.IGU#  YUPIGIl/  XCT.3P/  YE'nP 

ccMcnr,  /tchiz/  u.X/  my 

CfJMM'JM  /PLTCTP/  MAPID,  CCALZ/  MPLTIIT/  PLTHT(IO)*  ICDS,  IMP 
CUIIMOM  /r:g/  Hnv.':i:(A)00> 
rOMflOM  /AGMZTF/  n.lllf! 

COMKHM  /BOX/  POY.'IIM/  BCXIIAX/  BJYMZtl/  PO'^MAX/  ZDilX 
COr.llGM  /OTULF.;/  rPrrC.V/  ITCTAY/  V\'J\,  XMAX/  YillU/  YMAX/  NThPEZ/ 
''  “GF.I?/  IZAVZ 

CnMf'.OM  /VCGZMn/  jvcc”' 
riAL  l'!owh;(4ooo) 

CUMMGM  /ZZ/  2GLD(4JOO}/  ZMOIKACOO)/  2L, . TCP, ( 4000) 

Call  r.rMCHK 

PUT  r^'A  cllvaticm:  imtg  dip.cct  acc^c:  f:l»  (umit  s). 

Iir'^irr  r.-JUMPAR','  CflUF.DItlATC:  to  CORF.LCP.lM.r  UITM  PMA  CCORDDiaTL-, 

*t::'r.i:c  jprcr.vcF.  ..md  pcmp  gbc.  cl:v,\":oii  (hubs  mo^  cmcludco) 


CALL  :i:PGlT(M."ArH:/  TGP.ID»  OPsIPP,/  XIIM/  XMAX/  YMIM/  YMAX/ 
»  MX,  MY/  XGr.rOM.-  YJFCCMj  I'CUP/  YGTOr) 

CALI  ?:ru.G;T(XGO:/Yr)DG/  CPI'.TM.V) 

«  CLVUGCIXUnCi  YGI  C/Y.UF.ir.M/YJr  CGtl/ Of.CPR) 
l.'F.CTLCt/  9'’n)  GPTP.AM 


PLOT  PF.CLIMCMAF.Y  CMFU,  (AXCC/  OnGJP.VuP,  C GUP.D ItUTES/  ZTC). 

PLOT  LO:  :MFU  pop.  SCAMC  CLOCCC’!'  to  DDCCr.VCP. 

r/.v:  :  jld  im  lug. hi;  -  ucc  :?  as  ist  zgld  fur  scams  wcst  or  mbs. 

sa'.'l  zla'^ep  -  use  :t  as  ist  zglp  fgr  scams  east  gf  mbs. 

IF  SAVZMG  MAP  Oil  PISK/  PUT  PF.ZLIM*  ZMPr),  IM  RL'C  ‘tOCl-4G13  f.  SAVE 

::tgt  -  ZM,'!  (XGc.c-xur.:GM+,ooooi)/GL:L  r )  *  cniDr,  •••  xof.igm  +  gpiuF 
call  f ltpfe {x:i:m/  xmax/  ym:?’./  ymax) 

IF  (ISAV:  .EP,  1)  CALL  SAVPFZE Z'GTaV) 
c*,LL  rr.zTiMC'i-c'-/  lhouh:/  mgwmZ/  /ZCOLd 
C..LL  FL’LfKLH  )'.JM:/  ZcOLl-l/  Z'CTAV) 

CALL  PLTLIKHrjWU:/  ICGLl#  I'GTA'.') 

IF  (ISAVC  .LQ.  :)  Call  WFITMSO/L'.ini.Mini)/  MY/  ICMLl-i) 

A- 2 


:f  Cr/.'/C  .CQt  1)  CALL  WP.ITHSO#  HOHHICD#  NY#  ICOLl) 

PU  200  I  -  1#  NY 

LUnWHKI)  ■  ZOLD(I) 

200  CQMTIMUC 

c  **♦ 

(;  TILL  imiiH:  C  PLOT  FOK  LACH  LINO  EAST  OF  OBSERVEP. 

C  *** 

ICljL2  -  ICDLI  ♦  1 
ir  <:C[)L2  .CT.‘  NX)  CO  TO  450 
['u  400  :coL  •  rc  jLz#  hx 
rn  300  I  •  1#  MY 

ZHLDd)  -  ZLATlR(I) 

300  CUMTIMUE 

CALL  RCAtMS(2/  ZMOlMl)#  MY»  ICUL) 

CALL  OMELNdCOL#  HOWHI) 

CALL  PLTLLKHOWHI#  ICQL#  ITGTAV) 

IF  dSAVL  .no,  1)  CALL  IIRITHSO#  HOWHKD#  NY#  ICDL) 

4C0  COMTIHIJE 
450  CUtlTIHUL 

C  *** 

C  FILL  linuill  t  PLOT  FOR  EACH  LINE  WEST  OF  OBSERVER 
C 

rCLlLM?.  -  *C0L1  -  2 
IF  dCnLM2  ,iT,i  1)  GO  TO  f  00 
ICOL  -  ICOLl  -  1 
DO  50C  ;  -  1#  MY 

ZLATEP.d)  -  LMOHHId) 

500  CnMTIMlIE 

DU  700  I  ■  1#  IC0LM2 
ICOL  -  ICOL  -  1 
no  AGO  J  -  :#  MY 

ZOLIXJ)  ■  ZLATCR(J) 

6Cn  CONTINUE 

CALL  RLAPHSiZ#  ZMUUd)#  NY#  ICOL) 

CALL  ONELNdCOL#  MQUHI) 

CALL  PLTLNCHOWHI#  ICOL#  ITGTAV) 

IF  tlSAVE  .£Q.  1)  CALL  URITHSd#  HOl/MKl)#  NY#  ICOL) 

700  CJN^IMUL 

800  CONTINUE 

'U  dcox  .NC,  0)  CALL  PRFRSN 
CALL  PLTPGE 
STOP 
C 

c  format: 

c  **♦ 

950  ror.MArc  llevatimn  or  terrain  at  obccrver  position  ■•#  fio.d 

CML' 
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SUBROUTINE  DCDHAX  (  X#  Y»  HREC) 


*** 

44i* 


c 


BY  HONTC  COIENAN 

RETURNS  X  AND  Y  COORDINATES  FROM 


DATA  RECORD  ALUNG  WITH  RECnRO  NO 


COMHON  /OMACOH/  lBUr(700)#  ITEMP(flO)#  ITEM 
DATA  XSCALE/  ItOO  /»  YCCALE  /  1,00  / 


CALL  UNPACK  (IBUF#  ITEMP#  10) 
NRCC-ICVT(lTCMP(A)i  3) 
X-FL0AT(ICVT(1TCMP<7)#2) )*XSCALC 
Y-FL0AT(ICVT<ITLMP(0),2) )*YCCALE 
RETURN 
END 


I 


ooonortrio 


SUBROUTINE  OCTID  (ISH#  ISER#  lED#  XCOR#  YCUR) 

*** 

BY  HONTE  COLEMAN 

DECODES  A  TAPE  ID  BLOCK 

ISH  -  MAP  SHEET  IDENTIFICATIQM 

ISLR  -  MAP  SERIES  IDEIlTinCATION 

lED  -  HAP  EDITION  IDENTIFICATION 

CDMIERS  IN  QRDCr.  S\l#  NU»  NE»  SE 

CnNHON  /DMACOH/  IBUr<700)#  ITENPC80)/  ITEM 
DINENSION  XCaR(4),  YCDR<4)i  ISHt2) 

CALL  UNPACK  (IBUF#  ITEHP#  7C) 

CALL  roTATP.  <ITENP<7)/  12) 

CALL  PACK  (ITEMP(7)#  ISHi  12) 
ISi:R-ICYT(ITEMP(6-»)/  6) 
irD-ICVT(lTCI1P(73)#  6) 
lYC-l? 

IKC»2r 

no  100  I-  ifA 

XCGr.(I)-rLnAT<ICVT<lTi-Hp{;KC),6))*XSCALC 
YCnP.(I)-rHUT(ICVT(:TENr{lvc),6))^YSCALC 
IXC-IXC4-12 
IYC-IYC4-12 
100  CONTINUE 
RETURN 

DATA  XSCALC  /  1,00  /*  YSCALE  /  1.00  / 

END 
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FUMCTION  lILVUDrCXllRC#  Y^lQSi  XllP.tGll>  YDP.ICh,  CF  IDP  ) 

4i4i4< 

BY  BARBARA  FJ.JL1MS 

riMD?  7HC  GBscrA'-r.  i'LLV/“:n:i 

*** 


(>:?i  Y2i: 


:3) 


(Xfis:,  Yi,  nnuK) 


(;7»>  YA,7'») 


C 


C 


C 

c 


c 


riMLtii  :Gfi  i'.TrcfK  A(’uu) 

PAmL  ITTSCFKAOOO) 

ctiiiinii  /sc/.ii:/  iix<  ii'' 

ILCRT  «  II!T<(XOB3  -  Xjr.IOU  4-  *JO')Oi)  /  GRinP.)  +  1 
IfACMT  «  :LlF7  +  1 

iFALdi;  -  zvr^{^•^r>,:  -  vunrcii  +  .oono.)  ;  oi.iPk)  4  1 

lABnvc  ■  UlL'h;  +  1 

n 
y,z 
>'.3 

XA 

Y1 
Yc 
Y3 
YA 

4141* 

CALL  r::Annr(2,  lftcciud#  iiy,  il:ft) 

CALL  F"Mnfi:(2/  FTCCMd),  tl’,'/  IRIOflT) 

4<4i4‘ 

zi  »  LFTrcncr-LLtiu) 

Zl  -  r.l  -  TVCGC.l) 

Z2  -  I TTCCI’dABllVi:) 

72  -  r.2  -  TWGCr?.) 

Z3  -  F.i^CCnCAG'lV!:) 

Z3  -  23  -  TVL0(7.3) 

ZA  -  r.TrcuCB'jLtr.j) 

•  ZA  -  TVLGC'A) 

4i4<4< 


xcir.ioti 

+ 

CLLFT  -  1)  - 

'•  Of.  ICR 

Xi 

X'lr.IGi! 

+ 

C RIGHT  -  :) 

•<  GriPFs 

XB 

YJr.lGt! 

•f 

CBuLDU  -  1) 

«  or,: DR 

Ynr.ICtl 

+ 

CAt  n’':  -  1) 

*  or.: DR 

Y?. 

y: 

OLiL'n  -  Zl  +  HY.ur:  -  xi)/r.f,’rF.)  *  (za- 
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5 


UP  ■  Z2  {tX0BS-X2)/GRlDRl  •  (Z3-22) 

ELVOBS  -  DDUN  ♦  {(YQBS-YD/GRIDR)  ♦  <L1P-0QWNJ 

RETURN 

END 


r»o  o  o 


CUDK.OUTIME  rPTATR  CA/M) 

BY  HQM*n  COLCIIAM 

cumvzrt:  ricLP  data  charactep.^  to  msplay  choc 

*** 

r:ncMs;;on  ib{64) 

DO  IJD  I-1»H 
j-iA{n*i 

100  COilTIflUC 
RLTURU 

C  FIELD  DATA  TO  DltlPLAY  CUDL  TABLL 

data  rc  /  0#  41,  ?o»  y»  0/ 


* 

3»  4# 

r>i  6, 

7# 

3# 

*’# 

lu# 

* 

Hi  12i 

13#  14# 

If# 

16, 

■<  T 
>  1  / 

lit 

It: 

If'i  20i- 

21#  r.  2 

#  23# 

24# 

2'' 
t*  • 

t  26# 

*> 

42#  36#  3-,  5 

D#  44# 

# 

■)#  4  3# 

* 

3f>#  4 

1#  ?1# 

0#  , 

0# 

0# 

4b#  0# 

>«< 

27# 

26#  29# 

30# 

31# 

32# 

33#  34# 

Ik 

3r# 

36#  ro#  c# 

49# 

47#  C#  0 

or>v“>  no<^  ooooo 


SUDF.nUTIME  rRCSCHUTGT#  YT0T>  IFCSYH) 

■kill* 

FY  BARBARA  BRUQHE 

GETS  nUMBRATOR  AND  OEUOHINATOR  FOR  CALCULATING  FkACTIUN  UC  ^REA 
nr  IMTCREST  SEEN  BY  OBSERVER  -  CONSIPERIHG  VARIOUS  TARGET  HEIGHTS 

COMMON  /PLTSTF/  MAPIO#  SCALE#  MPLTHT#  PLTHTUO),  lOBS#  INK. 

CUMMON  /ASNSTF/  RNUM(lO)#  OCUnM 

CUlUinU  /BOX/  BOXMIN#  BOXMAX#  BOYHIN»  BOYMAX#  I  BOX 

IF  YOU  AREN'T  III  AREA  OF  INTEREST#  DON'T  DU  ANYTHING  -  JUST  PlTURM 

IF  ((XTGT  .LT,,  BOXMIN)  ,UR,  (XTCT  .GT.  DQXMAX)  .UP.* 

*  (YTGT  *LT..  BtlYHIIl)  .OR.  (YTGT  .CT.  BUYMAX))  RETURN 

4iik4i 

:r  YOU  CAN'T  SEE  ANY  TARGETS#  ADD  1  TO  DLNnfllllATUR  THEM  RETURN 

:r  cfgsvm  ,gt.  fiPLTUT)  denom  ■  dehum  +  i 
ir  (irCSYM  .GT,  HPLTliT)  RETURN 

INCREASE  MUliERATOP.  FUR  CAC!!  TARGET  IICICIIT  THAT  CAM  BF  St  EM* 
INCREASE  OLNOIIIIIATUR#  THEN  RETURN 

***■ 

PLl  100  r  -  IFCSYH,  NPLTHT 
RNUMd)  ■  RNUH(I)  +  1 
IOC  CGNTIMlir 

ntnuM  •  DlMun  +  i 

RETURN 

ullD 
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c  4ii|** 

C  PY  BARBARA  BRnOME 

C  CALCULATEr  HEIWIII  FOP.  OCAUr  CLPCEST  TU  ODrtP.VcP. 

C 

cnrMnn  rizt  znLr'(Aooo)#  2Mge'(4ooo)j  zLATi:r,(4ooo) 

COMH.Jn  /AREA/  XflRlCIi#  YtiRIClU  XETGP/  YSTOP 

rriiinofi  /see/  xobs#  ygbs/  mubs#  obtp.aii 

CnHHOM  /SCAIlS/nXi  HY 
CnUMdM  GP.rDF, 

riMEiirruM  nawnjd) 

REAL  LHUllHKl) 

C  ♦** 

C  FIUP  SCAtlS  CLilSE'T  TO  nBSEP.Vir.*  LEFT  IS  HUD/  RIGHT  11  ZMOl! 

c  *** 

ICULl  -  IMT(<XnDS  -  XnniGM  ♦  ♦OOUOl)  /  GP.IDF.)  +  ?. 

GALL  RLAnMr(n,  zfinud),  hy/  icnti) 

CALL  RlADMS{2/  ZILPd)/  tlY,  ECOLl-l) 

c 

c  Finn  HOijHi  enr.  pts.  ♦.bovt  oL'Ci-P’.’cr.  nti  bhth  scams 

c  *** 

C  ♦-*<*  FIRST  TWO 

ITOTl  *  ItlT((YUP.S  -  Y']R:CM  ♦  ,00001)  /  CniOR)  +  2 

:r  (iTGTi  »cT,  n“)  cu  T'j  izr- 

XTGT  -  xar.icN  +  (iCHLi  -  1)  ♦  nr.iDi,  . 

YTOT  -  YORICM  +  (ITOTi  -  1)  *  GRIOR 
HUWtlldTCTl)  -  HIT  (TVr,C(  ZHOU <  ITCTl) )  ) 

ZLATEF.dTGTl)  -  IHT(  ZtlDlH  rTGTl) ) 

LflUllHIdTGTl)  -  IllTdVECdJLndTGTl)  )) 

ZOLDdlCTl)  -  IHKZOLDdTCTl)) 

ITGTPl  •  ITCTl  +  1 
IF  CTCTPl  .GT,.  MY)  CO  TO  12S 
XTCT  ■  XTCT  -  CP.IPR 
C  ***■  ALL  OTflEP.C 

DU  IJO  :TGT  -  ITGTPx/  M'' 

YTGT  -  YTOT  +  GP.IDR 
C  ***  RIGHT  IF  OBrrRVCR 

XTGT  -  XTGT  +  CRIOf. 

ZP  •  ZPRIIIKXTO':’/  YTGT/  ITGT/  ZULLi  TLATtF.) 

Iiainil  CTCT)  -  ZP  -  ZtlOUCTGT)  +  TVECdMUlJCTOT) ) 

ETV'LO  -  ir;T(TVrG{Zt!UV.'(ITCT))) 
unUfllCTGT)  -  AMAXKlIOllllIdTGT)/  CTVEC) 

EZnnw  -  lMT(ZllUH{iTCT)) 

ZLATCP(ITCT)  «  AMAXKZP/  EZtlUW) 

c  +♦*  LIFT  nr  oBrcp.VLP. 

XTGT  -  XTCT  -  CRIDf. 

ZP  •  2PRlME(XTf,T/  YTGT/  ITGT,  ZLATER/  ZOLD) 

LIIJUMKITCT)  -  ZP  -  (ZnLDCTGT)  -  TVECt  ZUlp  ( I  TOT) ) ) 

LT'TG  -  :M7(T''rc{znLr'<ircT) ) ) 

LIIUL-HKITCT)  ■  AdAXiaHOL'HI  ( ITGT)  /  lTVlG) 

LZOLD  -  iriTCZOLOdTCT) ) 

ZOLC(ITCT)  -  AMAXKZP,  LZULD) 

100  COtlTIMUE 
C  ♦♦♦ 

C  riHD  HCIUHI  FOR  {-TZ*  PELLIU  OR  LE'UL  l.'ITlI  UfS,  UM  BOTH  SCANS 

C  ♦♦♦ 
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.:i  '  .-V..  ■  A-W' 


C  f-IF.:5T  Tl/U 

IZl'  ITCTMl  ■  ITCTl  -  1 

XTCT  ■  XOr.ICM  *  CCOLl  -  1)  *  GRIOK 
YTGT  -  YORIGM  ♦  (ITGTMl  -  1)  ♦  CRIDP. 

C  PUIMTi  HORIZOHTAL  TH  UBS*#  IF  THEY  EXIST 

IF  (ABS(YTCT  -  Y3BS)  .GT*  *00001)  GO  TO  1?0 
lini.'llICTCTIll)  -  II1T(TVLG(ZIIQU(ITGTM1))) 
x^LATLr.CTCTHl)  ■  ItlT (ZflHlK  ITGTMl ) ) 

LIIOUHrCTGTI'.l)  -  IMKTVCGdOLndTGTHl)  ) ) 

ZULD(ITGTMl)  -  IIlKZJLDdTr.THl) ) 

ITGTMl  -  ITGTM:  -  1 

XTGT  «  XnPlCM  +  (ICHLl  -  1)  *  GPIDP. 

YTGT  -  YORICM  ♦  dTGTfll  -  1)  *  CRIDP. 

If  (ITGTHl  «LT,i  1)  GU  TD  300 
C  *-'♦  nr.GT  TWO  Bt'LUW 

list  HUllHIdTGTni)  -  :flT<TV’CCd;iOU(ITGTfil}) ) 

Zl.ATCP.dTCTni)  -  IMT(ZliLlWdTCTMl)) 

IHaUlirdTGTIll)  -  :MT(TVCC(ZOLndTGTni)  )) 

ZOLLdTCTIll)  -  ::iT(ZGLDdT6T!!l)  J 
ITCTMP  -  ITCTHl  -  1 
:TGT  ■  ITGTtli 
XTGT  -  XTGT  -  CF IDR 
IF  (ir^THP  *lT,i  1)  CG  TU  300 
C  *«»>*«  /,Ll.  OTMCF.: 

DU  ZOO  J  «  1»  :tgtii2 
ITCT  -  ITCT  -  1 
YTCT  -  YTGT  -  CF.IDR 
L  f.icuT  or  udrj'.vcp 

XTCT  -  XTGT  +  GF.IDI. 

ZP  ■  ZPF.IHEFXTGT#  YTGT^  ITCT#  ZOLD,  ZLATER) 
llUWIIIdTCT)  «  ZP  -  (ZtinudTCT)  -  TVtO(ZIIOW(ITOT))) 
LTl'CG  -  :MT(T”t:C{Zliai.'dTGT)  ) } 

MlHIIIIdTCT)  -  AMAXKlinUHIdTCT)#  LTVCG) 

EZMnw  •  IHT(ZlinwdTOT)) 

ZL/.TLPdTCT)  -  AMAXKZP#  uZMUU) 
r  ILFT  OF  nFirCPVLP 

XTCT  ■  XTCT  -  GRIDR 

zr-  »  ZmilFO'^GTi  '"TGT*  I^GTi  ZLATLR/  ZOLD) 
LMD'.JIlIdTC")  =>  Zf  -  (ZOLrdTGT)  TVLG(  ZULDdTGT)  ) ) 
•ITVuG  -  1I1T(T''FC(Z'1LD(:TGT))) 

I  ilUWMKITGT)  »  AflAXKLinWIirdTCT)/  ZTVLC) 

LZ01.D  -  :fr(zuLn(;-'OT)) 

ZOLl'CTCT)  -  AMAXKZ.P/  CZOLD) 

CClMTIMlir 


SUBROUTINE  GQ2IN2(RLITL#BIG#LITNAH#BIGNAM) 


C  *** 

C  BY  BARBARA  BRQOHE 

C  HAKES  SURE  RLITL  GOES  INTO  BIG  (WITH  NO  REMAINDER), 

C  LITNAM  ■  *RLITL‘j  BI6NAM  -  'BIG'. 

C  *♦* 

IF  ( (AINT{BI6/RLITL)*RLITL)  .NE*  BIG)  WR ITE (6, 1^0 ) 

♦  BI6N AH/ LITNAM# BIGNAH^BIGi LITNAM/ RLITL 
IF  ( ( AINT(BIG/RLITL)*RLITL)  .NE*  BIG)  STOP 

RETURN 

100  FORHATdX/  A6/  •  SHOULD  BE  DIVISIBLE  BY  •/  A6/  // 

♦  IX#  A6/  •  -  »/  F10.3/  // 

♦  IX/  A6/  '  *  •/  FI0.3) 

END 


o  o  o  r>  r> 


FUNCTION 


ICVT  (IBUF#  N5 


♦♦♦ 


BY  MONTE  COLEMAN 

RETURNS  INTERNAL  INTEGER  VALUE  ASSOCIATED  WITH  M  6-aiT  BYTr 
STORED  ONE  PER  WORD  IN  ARRAY  IBUF.  *• 


DIMENSION  IBUF(l) 

INTEGER  SHIFT 

IF  (N  .GTt  10)  GO  TO  200 

ITEM-0 

DO  100  I-1>N 

ITEM-SHIFT! IT£H>6).QR,IBUF{I) 

100  CONTINUE 
ICVT-ITEM 
RETURN 

200  PRINT  210#  N 
STOP 

210  FORMAT  (•  ICVT/ERROR . LARGE  N#  N-'#  110) 

END 


o  o  o  o  o 


SUBROUTINE  XHNHX(  IVALUE#  XNINi  IHAX#  VALNAH#  HINMAH#  MAXM  AM  ) 

*** 

BY  BARBARA  BRQQHE 

HAKES  SURE  XNXN  IS  LESS  THAN  XVALUE  IS  LESS  THAN  IMAX, 
VALNAH  ■  ‘XVALUE*#  HXNNAH  ■  »XMXN'#  MAXNAM  -  ‘IMAX*. 

4>*4> 

IF  (XVALUE  .LE.  IHIN)  WRITE(6#1L0)  VALNAH#HINNAM# VALNAH# 

*  IVALUC/MINMAM/IHIN 
IF  (XVALUE  .GE.  XHAX)  WRXTE(6#200)  VALNAH# MAXNAM# VALNAM# 

*  IVALUE#MAXHAM#IMAX 
IF  (XVALUE  .LE*  IMIN)  STOP 

IF  (XVALUE  .GE*  XHAX)  STOP 
RETURN 

lOU  FORHATdX#  A20#<  SHOULD  BE  GREATER  THAN  •#  A2w/# 

*  IX#  A20#*  •  •#  XIO#  /# 

*  IX#  A20#*  -  •#  IIP) 

200  FORHATdX#  A20#*  SHOULD  BE  LESS  THAN  •#  A20#/# 

*  IX#  A20#>  •  •#  XIO#  /# 

IX#  A20#*  ■  •#  IIP) 

END 


A-X4 


r>  o  o  o  o 


FUNCTION  LFRACKZ/K) 


1 


2 


RY  ARTHUR  GROVES 

THIS  FUNCTION  RETURNS  TEN  TIMES  THE  FRACTIONAL  PART  OF  Tllif. 
K-TH  (OF  THREE)  ELEVATIONS  RACKED  INTO  THE  WORD  Z. 


A-AINT(Z) 

C-A 

IF(K.Ea.l)GQTO  I 

A«lOOww»*(Z-A) 

C-AINT(A1 

IF(K.Ea.2)G0T0  1 

C»10000*^(A-C) 

C-AINT(C) 

C-.I*C  ♦  .001 
A-C-AtHT(C) 
a— ,05 
DO  2  1-1^4 
B-B+.l 

IF(A.GT,B)GQTO  2 
LFRACT-I-1 
RETURN 
CONTINUE 


END 


A- 15 


':!:2is:ii:essxsz!:s3, 


o  r>  o  r»  o  o 


orir>r>  oooono  r>or>o 


4<*4< 


SUBROUTINE  HERGIT (NTAPES#  TGRID#  6RI0R/  XNIN#  XNAX# 
$  NY#  XORION#  Y0RZ6N#  XSTQP#  YSTCP) 


YHINp  YMAXiNX# 


«4i4> 


BY  ARTHUR  GROVES  AND  BARBARA  BROONE 

MERGE  UP  TO  12  DMA  TAPES#  ELEVATIONS  GO  IN  DIRECT  ACCESS  ARRAY. 

OIKENSION  HINDEXC^OPI) 

DIMENSION  BUF(40C0) 

DZHENSIUN  1(12)#  XQ(12)#  YO(ia}#  RfBOOOb  XC(4)#  YC(<>t 
DIMENSION  XL(IZ#M#YL(l2#5i#:4CN(4)#YCM(6)#C(A) 

CGMHrJH  /VEGIHD/  JVEG 

DATA  (C<IJ#:«l#4i/ •SOUTHWEST*# ‘MORTHWEST*# 'NORTHEAST*# •SOUTHEAST*/ 
CALL  OPENHSUi  MZNDEX#  4001#  0» 

COMPUTE  THINHIHG  CONSTANT 

ISKIP  •  AINT(GR:DR/TGRID  ♦  all 
IBZG  •  lOOOOOUOOO 

READ  HEADER  RECORDS.  COMPUTE  INDIVIDUAL  AND  COMBINED  TAPE  START 


IED#  XC#  YC) 


5J 


xOO 

#i4i#i 


XSTART  •  1000000000. 

YSTART  •  1000000000. 
i>Q  100  I-i#NTAPES 
K  «  I  ♦  6 

IP  (I  .GT.  6)  K  -  I  ♦ 

CALL  R0DATA(t</ ICQOE) 

CALL  RTYPRCdCQDu#  1) 

CALL  OCTIOCISH#  ISER# 

00  50  N«ltf4 

XL(I#N)-XCU:) 

YL(l#N»"yC(NI 
CONTINUE 
XLa#5J*XC<iJ 
YLCI#5)-YC(1) 

XQ(I)  •  AMINKXCd)#  XC(Z)) 

YO(I)  •  AHIN1(YC(X»#  YC(4)) 

XSTART  «  AMIMl (XSTART#  XOCIM 
YSTART  -  AMINKYSTART#  YOd)} 

CALL  ROOATA(K#  ICOOE) 

CALL  RTYPRCdCODE#  2) 

Ldl'  ■  0 
CONTINUE 

CHECK  TO  SEE  IF  EACH  CORNER  OF  THE  AREA  TO  BE  MERGED 
IS  ON  ONE  OF  THE  TAPES. 


*** 


XCH(l) 
XCHiZ) 
XCM(3) 
XCM(4) 
YCM(X) 
YCM(2> 
YCM(3) 
YCH(4) 
NCQT-0 
DO  179 


•  XMIN 
-XMIN 
-XHAX 
-XMAX 
•YHXN 
-YHAX 
-YMAX 
•YMIH 

K*l#4 
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t  d 


1  '4C 


O  O  o  o 


OQ  ISO  I«1«NTAPES 
N(W«0 

OQ  125  N«l«4 

0«XCM(K)«(YL(XjN)-YL(Z»N^I)) 

$  -YCM(K)*CXm#M)-XL«I»Nn)» 

$  -XL{I#Hn)PYL(I#N)^XL(IiH)*YLC  I»N+l) 

IF(D,LE«0«)NOR«HOR4>I 
12f  CONTINUE 

IFiNOR«EQ«4)GQTO  175 
15U  CONTINUE 

NCOT-MCOT41 
URITE(6«A0)C(K) 

175  CONTINUE 

C  #4"# 

C  CONFUTE' Y  COORDINATE  OF  1ST  POINT  ON  EACH  KERGED  SCAN  LINt 
C 


Y  ■  YSTART  -  TGRID 
DO  20C  I-1»IBIG 

IF  ((Y  •LE.  YMINJ  oAND.  (Y+TGRID  ,GT.  YHIN))  GO  TO  300 

Y  •  Y  ♦  TGRID 

2yO  CONTINUE 

30 t  YQRIGH  -  Y 
I  SAVE  ■  I  -  1 
C 

C  CGHPUTE  fSUHBER  OF  POINTS  ON  EACH  MERGED  SCAN  LINE 

C  ♦♦♦ 

DO  I-2#IBIG 

Y  •  Y  ♦  GP.IDR 

IF  ((Y-GRIOP  YMAX)  .AND*  (Y  .GT.  YHAX))  60  TO  56;> 

400  CONTINUE 

500  NY  -  I 
C 

C  SET  UP  LOOP  TO  CYCLE  THROUGH  INCREASING  VALUES  OF  X 
C  *** 


*** 


600 


X  ■  XSTART  -  TGRID 
XLAST  •  XSTART  -  GRIDR 
NX  •  j 

DO  900  I«1^IBIG 
JIM  ■  0 
X  •  X  ♦  TGRID 

IF  (X  .GT.  XHAX4-GP.I0R)  GO  TO  IGuO 

CYCLE  THRU  THE  TAPES  PICKING  OFF  DATA  FROH  EACH  TAPE 
THAT  CORRESPONDS  TO  THE  CURRENT  VALUE  OF  X 

DO  600  J^1#NTAPES 
K  ■  J  ♦  6 

IF  U  ,GT,  6>  K  •  J  ♦  7 

IF  ((X  ,LT.  XQ(JI-.l)  .OR.  (L(J).EQ.in  GO  TO  6^3 
CALL  RODATAIK#  ICODc) 

IF  (dCOOE.NE.’EOF' I  .AND.  ( ICOOE.NE,  *  cOI  •  J  )  GQ  10  6 
L(J)  -  1 
GO  TO  800 

IF  (ICCjDE  .LO.  ’DATA')  GO  TO  700 
lUNIV  -  J  ♦  6 
WRITE C6/ 10)  lUNIT 
STOP 
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■  v--' 


■■frT  I'V' 


700  IF  ((X.LT.XHIN-TGRID)-nR.(X-XLAST.LT.GPIDR-i.))GU  in  , 

JIH  ■  1 

CALL  OCDHAXiXX#  YY>  NREC) 
lY  ■  YY 

M  ■  AINT{{YQ(JJ-YSTAP.T)/TGRID  ♦  .1) 

CALL  STOATAtR^H't’IYnieOOOiJVEG) 

00  L  CONTINUE 

C  *** 

C  THIN  AND  RECORD  THE  ELEVATION  DATA  FOR  THIS  SCAN  LINE. 

C  COMPUTE  NUMBER  OF  SCANS  AND  X  ORIGIN  OF  MERGED  DATA. 

C  *** 

IF  (JIH  .EQ.  l>)  GO  TO  9o0 
MX  ■  MX  ♦  1 

IF  (NX  .EQ.  1)  XQRIGM  -  X 
XLAST  -  X 

MSAVE  •  ISAVE  ♦  (MY-1)*ISKIR 
ICQUNT  -1 

DO  850  J-ISAVEi  NSAVE/  ISKIP 
BUFdCOUNT)  ■  R(J) 

ICOUNT  ■  ICQUNT  ♦  1 
050  CONTINUE 

CALL  WRITMS(2/  BUF(l),  NY>  NX) 

9l'U  CONTINUE 

C 

C  WRITE  IDENTIFYING  DATA  FOR  MERGED  TAPE 
C 

XSTilP  ■  XQRIGN  +  GRrDR*FLaAT(NX-l) 

YSTOP  -  YQRIGN  4  GRIDR*FLaAT( NY-1 ) 


C 

c 

c 


WPITE(6/20)  MX>  MY/  XQRIGN/  YQRIGN/  XSTQP/  YSTOP/  GP.IDP 
RETURN 

*** 


FORMATS 

*♦* 

10  FQRMATC  SOMETHING  NON-STANDARD  ABOUT  RECORD  FORMAT  QN  UNIT*/  II- ) 
21  FaRMAT(///  •  MERGED  ARRAYS  // 


$ 

i 

$ 

1 

$ 

i 

$ 


NUMBER  OF  SCAN  LINES  -•/ 

NUMBER  OF  POINTS  PER  SCAN  LINE  •»/ 

X  COORDINATE  OF  FIRST  SCAN  LINE  •'/ 
Y  COORDINATE  OF  FIRST  POINT  ON  EACH 


X  COORDINATE  OF  LAST 
Y  COORDINATE  OF  LAST 
ARRAY  GRID  -•/ 

30  FORllATC///  '  RECORDS  110/  // 

S  (•  S  1DF7,*)) 

40  FORMAT! ZIX/AIO/ 'CORNER  OF  AREA 
S'  ) 


Ilv^//» 

I  Iv  » /  / 
FiO.l  / 

SCAN  LINE  -'/Flw.lW; 

F  iv.  .  i  /  / » 

POINT  ON  EACH  SCAN  LINE.  -•»  FlO.iW/ 

Flw.l) 


SCAN  LINE  -S 


NOT  ON  ANY  TAPE  -  RUN  ABORTED 


END 


-.I, 
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SUBROUTINE  ONELNdCQL#  HOWHX) 

BY  BARBARA  BRQQHE 

GIVEN  ARRAYS  ZOLO  (HAX  PREVIOUS  INTERRUPT)  E  ZNQW  (CURRENT  DMA 
ELEVATIONS  PLUS  VEGETATION  HEIGHT  AND  VEGtTATION  CODE) 
CALCULATES  FOR  ONE  SCAN  HOW  HIGH  TGT  MUST  Bb  RAISED  TO  RF 
SEEN  E  UPDATED  HAX  INTERRUPT  (ZLATER). 

COHHON  /ZS/  ZOLD(AOOO)#  ZNOU(AOOO)>  ZLAT0R(A^'L  ^ ) 

COHHON  /AREA/  XQRIGN#  YQRI6N#  XSTOP,  YSTOP 
COHHON  /SEE/  XQBS#  YQBS>  HOBS/  QBTRAM 
COHHON  /SCANS/  NX#  UY 
COHHON  GRIDR 
DIHENSION  HQWHKl) 

C  *4)# 

C  SPECIAL  CASEt  OBSERVER  ON  A  SCAN  LINE 

C  AAA 

XT6T  ■  X0RI6N  *  (ICQL  -  1)  A  GRIDR 

IF  (ABS(XTGT  -  ( X0BS-6RIDR) )  *GT.  .OuOul)  GO  TO  30 

ITGTl  •  INT((YQBS  -  YQRIGN  ♦  .OOuOl)  /  GRIDR)  A  2 

YT6T  -  YQRIGN  ♦  (ITGTl  -  1)  *  GRIDR 

ITGTHl  ■  ITGTl  -  1 

ITGTH2  ■  ITGTl  -  2 

IF  (ITGTl  .GT,  NY)  GO  TO  25 

HQWHKITGTl)  ■  INT( TVEG(2H0W( ITGTl )) ) 

ZLATERdTGTl)  •  INT(  ZNOW(  ITGTl) ) 

25  IF  (ITGTHl  .LT.  1)  GO  TO  50 

HQUHKITGTMI)  •  INT  ( TVE6(  ZNQW(  ITGTHl )) ) 

ZLATERdTGTHl)  -  INT(ZNOWdTGTHi) ) 

IF  dTGTH2  •LT.  1)  GO  TO  50 

IF  (ABSiYTGT  -  { YDBSAGRIDR) )  .GT.  .Ok^Oi^l)  GO  TO  Du 
HQWHKITGTHE)  •  INT(TVEG(  ZNOW(  ITGTH2 ) ) ) 

ZLATER(ITGTH2)  -  INT(ZN0WdTGTH2)) 

C  AAA 

C  POINTS  ABOVE  OBSERVER 

C  AAA 

50  ITGTl  •  INT((YOBS  -  YQRIGN  ♦  .00001)  /  GRIDR)  +  2 

IF  (ABS(XTGT  -  (XOBS-GRIDR) )  .LT.  .uuOo: )  ITGT*  -  ITGTo.  +  ^ 
XTGT  •  XQRIGN  +  (ICOL  -  1)  A  GRIDR 
YTGT  -  YQRIGN  +  (ITGTl  -  1)  A  GRIDR 
IF  (ITGTl  .GT.  NY)  GO  TO  150 
DO  100  ITGT  •  ITGTl#  NY 

ZP  •  ZPRIHE(XTGT/  YTGT#  ITGT#  ZOLD#  7.LATER) 

HOWHKITGT)  -  ZP  -  (ZHOW(ITGT)  -  TVEG(ZMOW(  ITGT) )  ) 

ETVEG  -  INT(TVEG(ZNOM(ITGT)) ) 

HOWHKITGT)  -  AHAXl  (HOWHKITGT) #  CTVEG) 

EZNQW  -  INT(ZNQW(ITGT) ) 

ZLATERdTGT)  -  AHAXKZP#  EZNOW) 

YTGT  •  YTGT  ♦  GRIDR 
ILO  CONTINUE 
C  AAA 

C  POINTS  BELOW  OBSERVER 

c  AAA 

150  ITGTHl  •  ITGTl  -  1 

IF  (ABS(XT6T  -  ( XOBS-GRIDR ) )  .LT.  .OtOOl)  ITGTHl  -  ITGTx  -  3 
YTGT  -  YORIGN  ♦  (ITGTHl  -  1)  A  GRIDR 
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IF  (ABS(YT6T  -  (YOBS-GRXDR) )  .LT.  .00001)  ITCTHl  -  IT6T1  - 
IF  (IT6TN1  .LT.  1)  GO  TO  300 
ITGT  ■  ITGTMl  ♦  X 

YTGT  •  YORIGN  ♦  tXTGT  -  1)  ♦  GRIOR 
DO  200  X  ■  1#  XTGTNl 
ITGT  •  XTGT  -  1 
YTGT  •  YTGT  -  GRIDR 

ZP  •  ZPRXHE(XTGT«  YTGT«  XTGT.  ZOLD.  ZLATER) 

HOUHKXTGT)  -  ZP  -  (ZNOUdTGT)  -  TVEG( ZNOU(  ITGT) ) ) 

ETVEG  •  INT(TVEG(ZNQW(ITGT))) 

HOWHX(XTGT)  •  AHAXKHOWHXdTGT) .  ETVEG) 

EZNOU  •  INT(ZNQU(ITGT)) 

ZLATER(ITGT)  -  AHAXKZP#  EZNQW) 

200  CONTXNUE 
300  RETURN 
END 
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SUBROUTINE  PLTICS < IXRY AX#  BOTRLF#  START#  ST  IP) 

*** 

BY  ARTHUR  GROVES  AND  BARBARA  nROOHE 
THIS  ROUTINE  PUTS  LARGER  TIC  HARKS  ALONG  AN  AXIS 
IXRYAX  -  INDICATES  WHETHER  HARKING  X-  OP  Y-AXIS  (i  -X  ,  ^-Y) 
BOTRLT  -  GIVES  THE  aOTTHH  OR  LEFT  (DE “EllDIMG  OS  IXRYaX)  T'< 
IF  IXRYAX-C#  GIVES  BUTTON;  ZlZ'l  GIVI  :  FT 

start  -  STARTING  C  IORDINATl  1.,  \X  IS 
STOP  -  STOPPING  CQQRCIt.AT;;  M,i  AXIS 

DIMENSION  EKS{2)#  WYC(2) 

C  *** 

C  TO  LABEL  X-AXIS 
C 

IF  (IXRYAX  ,EQ«  i)  GO  TO  2^. 

WYt(i)  •  BQTRLF 
WYE(2I  ■  BQTRLF  ♦  60. 
on  IOC  I  >  1#  1000 

EKS(l)  •  START  ♦  iOco.^FL OAT ( I-* ) 

CKS<2)  -  EKSCl) 

IF  ((EKS(l)-STOP)  .GT,  5u.)  F.tTUR*! 

CALL  PLTDTSd#  0#  £KS>  WYu#  a#  ■.  ) 


CONTINUE 

P.tTURN 

c 

c 

TO  LABLl  Y-AXIS 

c 

Z<j*j 

EKS(l)  •  BQTRLF 

:KS(2)  ■  BQTRLF  +  60, 

DU  300  I  •  1#  IwOO 

WYE(I)  •  START  ♦  10Lv,'.*FL  )AT  ( I-.) 
WYc(2)  ■  WYE<1) 

IF  ((WYE(l)-STQP)  ,GT,  50.)  hlTuni. 
CALL  PLTDTSd#  0#  EKS#  V/YC#  2/  '  ) 
COIITINUE 
Rf  TURN 
END 
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SUBROUTINE  PLTLNCHOWHI#  ICOL#  ITGTAV) 

C  ♦♦♦ 

C  BY  BARBARA  BROOME 

C  PLOTS  ONE  SCAN  LINE  OF  HaWHX*S  (ODD  COLUMNS  UP  OTHERS  DOWN) 

C  IF  KEEPING  UP  WITH  FRAC.  OF  BOX  SEEN/  CALL  NEC.  SUBROUTINE 
C 

DIMENSION  XPL0T12)/  YPLOTIZJ 
DIMENSION  SYHBOL(IO) 

DIMENSION  HOWHKACOu) 

COMMON  /ZS/  Z0LD(40U0)/  ZNOW(4000)/  ZLATcR(40uu) 

COMMON  /AREA/  XORIGN/  YORIGN/  XSTOP/  YSTOP 
COMMON  /SCANS/  MX#  NY 
COMMON  GRIDR 

COMMON  /PLTSTF/  MAPID#  SCALE/  NPLTHT/  PLTHT(IO)/  lOBS/  INK 
COMMON  /BOX/  BQXHIN/  BOXMAX/  BOYHIN/  BOYMAX/  IBOX 
DATA  SYMBOL  /•+>•/  */>•/  •♦>•/  »B>'/  'C>'/  •D>»/  '£>•/  •F> 

*  (->«/ 

SCLCNV  •  SCALE/100. 

HIGH  -  .ASAGRIDR 
WIDE  -  .ZBCBAGRIOR 
SIZE  •  .9PGRIDR  /  SCLCNV 
D  -  .sagriop 

IFIHODdCOL/  ^)  .EQ,  C)  D  •  -0 
XTGT  ■  XORIGN  +  { ICQL-1) PGRinR 
IF  (NPLTHT  .L£.  1)  GO  TO  AOO 
C 

C  CASE  li  MULTIPLE  TARGET  HEIGHTS 

C  GIVEN  SCAN  OF  HQWHI*S/  FIGURE  (FOR  EACH  Y)  WHICH  SYMBOL  TO  PLOT 
C  AND  PLOT  IT 

C  *** 

DO  300  I  •  1/  NY 

IF  (MQDdCQL/  2)  .EQ.  1)  J  •  I 

IF  (MQDdCQL/  2)  .EQ.  0)  J  •  NY  -  I  +  1 

YTGT  •  YQRIGN  ♦  (J-1)*GRI0R 

DO  100  IFGSYH  •  1/  NPLTHT 

IF  (ITGTAV  «EQ.  1)  HQWHKJ)  -  HQWHX(J)  - 

*  INT(TVEG(ZNOW(J) )) 

IF  (HOWHKJ)  .LE.  PLTHT (IFGS YH) )  GO  TO  Zifr 
lUU  CONTINUE 

IFGSYH  -  11 

200  IFdFGSYM  .GT.  1)  CALL  PLTSYH(SIZ£/  SYMBOL (  IFGS YH-1 ) /  i#./ 

*  XTGT -WIDE/  YTGT-HIGH) 

IF  (IBOX  .NE,  0)  CALL  FRCSOH(XTGT/  YTGT/  IFGSYH) 

300  CONTINUE 
RETURN 

C  *•¥•¥ 

C  CASE  21  YES/NO  HAP  (ONE  TARGET  HEIGHT) 

C  XPL0T(1)£(2)  AND  YPL0T(i)£(2)  MOLD  START/STOP  OF  O-O-V  SEGMENTS 
C  EVEN  COLUMNS  PLOT  DOWN/  ODD  COLUMNS  PLOT  UP 
C 

AOC  XPLOT(l)  -  XORIGN  +  ( IC0L-1)AGRIDR 
XPLaT(2)  ■  XPLOT(l) 

YPLOTd)  ■  0. 

YPL0T(2)  ■  0. 

DO  500  I-l/NY 


w.  "m|» 


I 


f 

ti 


IF  <M0D1ZC0L#2)  *£0*  0)  J  -  NY  -  X  4-  l 
YT6T  •  YONICN  ♦  {J-l»  *  GRIOR 
IF  (ITGTRV  tEQ.  X)  HOWHKJ)  •  HOWHZ(J)  - 

♦  XNT(TVEG(ZH(MtJn) 

IF  (HOWHZ(J)  «LE«  FlTHTCin  ZFGSYN  -  1 

IF  (HQUHZtJI  .GT,  PUTHTU))  IFGSYH  -  99‘)9<) 

IF  (IBQX  .HE.  0)  CALL  FRCSENCXTGTi  YTGT>  IFGSYH) 

IF  C(YPL0T(1)  «LE«  C«)  •AND*  (HOUHI(J)  .GT.  FLTHTd))) 

♦  YFLQTCD  -  YORXGN  ♦  (J-1)*GRIDR  -  D 

IF  (YPLQTCl)  «LE«  0.)  GO  TO  900 

IF  (HOWHZ(J)  tLE.  PLTHTd}) 

♦  YPLQT(2)  ■  YORIGN  ♦  (4-l)*GRIDR  -  0 

IF  (I  .EQ.  NY) 

♦  YPLQTC2)  •  YORIGN  ♦  (4-l)*GRIDR  4  D 

IF  (YPLQTI2)  •LE.  0.)  60  TO  900 

CALL  PLTOTSd#  C#  XPLOT#  YPLQT«  2,  0) 

YPLQT(l)  «  0« 

YPLQTiZ)  •  0. 

SOO  CONTINUE 
RETURN 
END 


! 


f 
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SUBROUTINE  PLTFRE(XHXN^  XHAXi  YHIH#  YHAX) 

BY  Barbara  broone 

PLOTS  preliminary  INPORMATION  -  BEFORE  LOS  CALCULATIONS 

4>4i4i 

DIMENSION  PKTS(2) 

DIMENSION  POBS(5) 

DIMENSION  PLB(6} 

DIMENSION  HID(2) 

DIMENSION  LABELPO)#  LABELI(4)«  SYMBOLdU) 

DIMENSION  LGN01(1)>  L6N02(3)/  LGND3(3)#  LGN04(2)/  LGH0^(2) 
COMMON  /PLTST'/  MAPID#  SCALE#  NPLTHT#  PLTHTdO#  lOBS#  INK 
COMMON  /SEE/  XOBS#  YOBS#  HOBS#  OBTRAN 
DATA  LABELP  /•BROOME  •#  •BLDG  392  •#  ‘XEAl?  •/ 

DATA  LABELI  /•BROOME  •#  •BLDG  392  *t  'XZAIT 

♦  •INK  PLEASE •/ 

DATA  L6N01  /•LEGEND  >•/ 

DATA  L6N02  /•POSITION  W*#  ‘HERE  TARGE*#  ‘T  >•/ 

DATA  LGND3  /'X  METERS  H'#  ‘IGH  OR  MOR*#  •£  >•/ 

DATA  LGN04  /•CAM  BE  SEE*#  •N#  WHERE  >*/ 


DATA  LGND5  /'SYMBOL  *# 

1 

*#  'X 

■  >•/ 

DATA  SYMBOL/'  >'#  '♦>'#  •/>* 

,  I4>« 

# 

*A>'# 

•B>'#  'C>*# 

♦  'F>*/ 

c 

c 

PLOT  AXES 

c 

4>4<4> 

PXMIN  ■  1000,  ♦  AINT(.OL'l 

4i 

(XMIN 

+ 

,Otl ) ) 

PXMAX  -  1000,  >•  AINT(.001 

4 

(XMAX 

» 

,001)) 

+  1000. 

PYMIN  •  1000,  *  AINT(,OOX 

4 

(YMIN 

♦ 

,001) ) 

PYMAX  •  1000,  *  AINT(.00l 

4 

(YHAX 

- 

•  ODD  ) 

+  looo. 

F  •  SCALE/100. 

XPA6E  -  (PXMAX  -  PXHIN)/F  *  32« 

IF  <XPA6E  .LT.  46.)  XPAGE  •  46, 

YPAGE  -  (PYMAX  -  PYMINJ/F  ♦  8. 

IF  « YPAGE  *LT.  21.)  YPAGE  •  21, 

CALL  RMNMX<XPAGE#45.99#3657,#«XPACE*#*TITLE  PLaT*#*12y  FT  MAX*) 
CALL  RMNMX(YPAGE#20.9#73,#  'YPAGE'#  'LABEL  SIZE'#  *29  IN.  MAX') 
CALL  PLTPGE 

IF  (INK  .EQ,  0)  CALL  PLTBEG(XPAGE#  YPAGE#  ,3937#  13#  LABELP) 

IF  (INK  .EQ.  1)  CALL  PLTB EG(XPAGE#  YPAGE#  ,3937#  13, LABELI) 

CALL  PLTSCA(6,#  6,#  PXHIN#  PYMIN#  F#  F) 

CALL  PLTAXSdOOC.#  1000,#  PXMIN#  PXMAX#  PYMIN#  PYMAX#  4) 

C  *** 

C  INCREASE  PLOT  TICMARK  SIZE  (TOP#  RT#  BOTTOM#  L) 

C 

CALL  PLTICS(0#  PYMAX  #  PXMIN#  PXMAX) 

CALL  PLTICSd#  PXMAX  #  PYMIN#  PYMAX) 

CALL  PLTIC5(0#  PYMSN-60.#  PXMIN#  PXMAX) 

CALL  PLTICSd#  PXHIN-6P.#  PYMIN#  PYMAX) 

C 

C  LABEL  AXES 

C  ♦♦♦ 

CALL  LABELAdOOp.#  1000.#  PXMIN#  PXMAX#  PYMIN#  PYMAX#  .001#  .OCl) 

C  ♦♦♦ 

C  PLOT  MAP  IDENTIFICATION 

C  *** 
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£NCaDE(20#  900/  HIDill)  HAPID 

CALL  PLTSYH(.3/  HID/  0*/  PXMIN/  PYHAX  ♦  .8*F) 

C  *** 

C  PLOT  aBSEPVER 

C  *** 

CALL  PLTOTSO/  0/  XOBS/  YDBS/  1/  0) 

C  *** 

C  PLOT  OBSERVER  COORDINATES  AND  HEIGHT 

C 

ENCODE (A2/  910/  PQBS(l))  XOBS/  YOBS 

CALL  PLTSYH<.3/  POBS/  0«/  PXHIN/  PYHIN  -  2.8*F) 

ENCODE  (29/  920/  POaSd))  HOBS 

CALL  PLTSYH(.3/  POBS/  0*/  PXHIN/  PYHIN  r  3.6*F) 

ENC0DE(29/  940/  POBSd))  IQBS 

CALL  PLTSYH(.3/  POBS/  0./  PXHIN/  PYHIN~4.4*F) 

C  *441 

C  PLOT  LEGEND 

C 

CALL  PLTSYM(.3/  LGNDl/  0*/  PXHAX+7.*F/  PYHIN+12.*F) 

CALL  PLTSYH<.3/  LGND2/  0./  PXHAX+7.*F/  PYHIN+li,2*F ) 

CALL  PLTSYH{.3/  L6ND3/  d»/  PXHAX47.*F/  PYHIN+10,4*F ) 

CALL  PLTSYH(.3/  L6N04/  0*/  PXHAX+7d4F/  PYHIN-f  9,6*F) 

CALL  PLTSYH(.3/  LGND5/  0./  PXHAX47.*F/  PYHIN+  8.e*F) 

DO  100  I  •  1,  NPLTHT 

CALL  PLTSYH(.3/  SYHBOL(I)/  0«/  PXKAX+8,4F/ 

♦  PYHIN*(8,8-I*,e)*F' 

ENC0DE(ll/930/PHTS(l)I  PLTHT(I) 

CALL  PLTSYH(.3/  PHTS/  0./  PXHAX+13.4F/ 

♦  PYHIN*(8.8-I4.8)*F) 

100  CONTINUE 

ENCODE (53/980/PLB(l) JPLTHT (NPLTHT J 

CALL  PLTSYH(.3/  PLB/  0./  PXHAX+8.4F/  PYHIN+( 8. 8-I*.8 )*F ) 

RETURN 

C  ♦** 

C  FQRHATS 
C  *4>4 

900  FaP,HAT(»HAP  ID  -  •/  AlO/  •>•) 

910  5  ORHATC  OBSERVER  C »/  'OORDlNATcS  */  *  (•/  F7,0/  •  /  •/  F8.L/  •)> 
920  FaRHAT( ’OBSERVER  H'/  ’EIGHT  •  ’/  FiO*2/  ’>») 

930  F0RHAT{F10.2/  •>•) 

940  FQRHATC ’OBSERVER  I'/  ’0  •  •/  110/  '>•) 

980  FORHAT(»~  NO  TARGET  <-’/F10.2/*  HETERS  HIGH  CAN  BE  SEEN>’) 
END 
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SUBROUTINE  PRFRSN 

♦♦♦ 

BY  BARBARA  BROOME 

PRINTS  THE  AVERAGE  AREA  SEEN  FOR  EACH  TARGET  HEIGHT 
COMMON  /ASNSTF/  RNUMdO)#  DENOM 

COMMON  /PLTSTF/  MAPID#SCAL£>  NPLTHT/  PLTHT(lit)/  lOBS/  INK 
COMMON  /BOX/  BQXMIN/  BQXMAX#  BOYMIN/  EOYHAX#  IBOX 
WRITE(6«  900)  BOXHIN#  BQYMIN/  BOXHAX/  BOYHAX 
WRITE(6»  910) 

00  100  I  -  l<  NPLTHT 

CFRACT  -  RNUH(I)  /  OENOM 
WRITE(6/  920)  CFRACTi  PLTHTCI) 

100  CONTINUE 
RETURN 

C  *** 

C  FORMATS 
C 

900  FQRHAT{///i  »  AREA  OF  INTEREST  DEFINED  BY  •/  /, 

♦  (•  (S  FlO.l*  *  *•,  FlO.l#  •)  •)) 

910  FORMAT{/i  •  FRACTION  OF  AREA  TARGETS  /i 

*  •  VISIBLE  TO  OBSERVER  HEIGHTS) 

920  FORHAT(/>  «  S  FlO.Di  7X#  FlO.l) 

END 
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SUBROUTINE  RDDATA ( IU> ICQDE) 


BY  MONTE  COLEMAN 

READS  ONE  BLOCK  FROM  A  DMA  TAPE  ON  UNIT  NUMBER  lU  AND  RETURNS 
ICQDE/  WHICH  INDICATES  THE  TYPE  OF  RECORD  READ. 

*** 

COMMON  /DMACQM/  IBUF (700)/lT£MP« 80)/ITEM 
INTEGER  TID/FID/DTA/EOF/EOI 

DATA  MSKA  /  777777777777000CCD00B//HSKB/  777777Uia  i  P,/# 

*  TID  /  00000000002000UODOOGB//FID  /  UUOCuOu  Jcl2L  wi  wU  .)  B/ , 

*  DTA  /  AlOLUctDCoUCUu&gOOOOB// EOF  /  0125A67.6l5uOJOOO(i juc 8 / 

*  EOI  /  0125513115DCcl-CDCDDLB/ 

BUFFER  IN  (lU/l) (IBUF{i)/IBUFI70u) ) 

IE0F«UNIT(IU) 

IFdEOF.GT.-DGOTO  2u0 

100  ITEH-IBUF(l) .AND. MSKA 
ICQDE* »  • 

IF<ITEH.EQ«TID)ICauE-‘TAPE  ID* 

IF(ITEM.EQ.FIO)ICODE-»FILn  ID* 

IF{ITEH.EQ.EaF)ICDDE-»EnF* 

IFdTEM.EQ.EOI  )  I  CODE -•EOI* 

ITEM-IBOFd),ANO.HSKB 

IF(ITEM.EQ.DTA)ICQD£-‘DATA* 

IFdCODE.EO*  '  •  )CaTO  ZAO 

RETURN 

200  IFdEQF.GT.O)GQTD  220 
PRINT  210/ lU 
ICaDE**EOF* 

RETURN 

210  FaRMAT<*  . . EOF  UH  UMIT'/I5) 

220  PRINT  230/ lU 
GOTO  100 

230  FORMAT!*  RDDATA/WARKING... ..PARITY  ERROR  ONUNIT'/IS) 

2A0  PRINT  250/ ITEM 
PRINT  260/ lU 

260  FORMAT!*  ERROR  ON  TAPE*/  110) 

STOP 

250  FORMAT!*  RDDATA/ ERROR.  ....  ILLEGAL  INDICATOR* /i^X/ » ICODE-  /»/02^’/ 

*  '/*) 

END 
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SUBROUTINE  RDNCHK 

*** 

BY  BARBARA  BROOME 

READ#  WRITE#  AND  CHECK  INPUT  FOR  SEEFAR  PROGRAM 

MAPIO#  SCALE#  INK  -  TELLS  WHAT  HAP  TITLE  TO  PUT  ON  PLOT, 

WHICH  SCALE  TO  USE  L  WHETHtR  TO  USE 
BALLPOINT  PEN  OR  INDIA  INK 

lOBS#  XOBS#  YOBS#  HOBS#  IDNTHV  -  OBSERVER  IDENTIFICATIUM# 

UTM  COORDINATES#  HEIGHT#  AND  INPICATrOM 
OF  WHETHER  TO  ATTEMPT  TH  EXPOSE  OBSERV.J 
NPLTHT#  IT6TAV  -  NO.  TGT,  HTS.  CONSIDERED  (LESS  OP  •  lo) 

AND  INDICATOR  FOR  WHETHER  HTGT  IS  HT 
ABOVE  V£G  (IF  1)  OR  ABOVE  LAND  (IF  u) 
PLTHTdl  -  I-TH  TARGET  HEIGHT  (METERS) 

XHIN#  XMAX#  YHIN#  YHAX#  GRIDR#  JVEG 

UTM  COORDINATES  OF  BOUNDARIES 
OF  AREA  TO  BE  CONSIDERED  £  GRID  TO  USE- 
INDICATOR  TO  ADD  VEGETATION 
VEG  e  URBAN»0#URBAN  ONLY-1# NOTHING-Z 
NTAPES#  TGRID  -  NUMBER  OF  DMA  TAPES  RcOUIRED  £  TAPE  GUD 
IBQX#  ISAVE  -  CALCULATE  VISIBLE  FRACTION  OF  BOX 

(1  -  YES#  0  -  NO) 

INDICATOR  FOR  WHETHER  TO  SAVc  HOWHI'S  ON 
DISK  FOR  FUTURE  COMPOSITE  CALCULATIONS 
(1  •  YES#  L  -  NO) 

HQWHI  AFFECTED  BY  ITGTAV  BtFORE  STORING 
BQXMIN#  BOXMAX#  BOYHIN#  BOYHAX  -  IF  IBOX  •  1#  DESCRIBES 

AREA  OF  INTcRtST 


*** 


COMMON  /VEGIND/  JVEG 
COMMON  GRIDR 

COMMON  /SEE/  XOBS#  YOBS#  HOBS#  OBTRAN 

COMMON  /PLTSTF/  MAPID#  SCALE#  NPLTHT#  PLTHT(ID)#  lOBS#  INK 
COMMON  /BOX/  BQXMIN#  BOXMAX#  BOYMIN#  BOYHAX#  IBOX 
COMMON  /OTHERI/  IDNTHV#  ITGTAV#  XHIN#  XllAX#  YHIN#  YHAX#  NTAPES# 
♦  TGRID#  ISAVE 

C 

C  READ  AND  WRITE  INPUT 

C 

READ  (5#9UP)  MAPIO#  SCALE#  INK 

WRlTn(6#9AC)  MAPID#  SCALE#  INK 

READ  (5#910)  lOBS#  XOBS#  YOBS#  HOBS#  IDNTHV 

WRITE(6#910)  lOBS#  XOBS#  YOBS#  HOBS#  IDNTHV 

READ  (5#96p)  NPLTHT#  ITGTAV 

WRITE(6#960)  NPLTHT#  ITGTAV 

READ  (5#920)  (PLTHT(I)#  I  ■  1#  NPLTHT) 

WRITE(6#920)  (PLTHT(I)#  I  -  1#  NPLTHT) 

READ(5#93L)XHIN# XMAX# YHIN# YHAX# GRIDR# JVEG 

WRIT£(6#930)XMIN#Xi’1AX#YMIH#YHAX#GRIDR#JVFG 

READ  (5#910)  NTAPES#  TGRID 

WRITE(6#910)  NTAPES#  TGRID 

READ  (5# 960)  IBOX#  ISAVE 

WRITE(6#960)  IBOX#  ISAVE 

IF  (IBOX  .EQ.  1)  READ  (5#92U)  BQXMIN#  BOXMAX#  BOYMIN#  BOYHAX 
IF  (IBOX  .EQ.  1)  WRIT£(6#92C)  BOXHIN#  BOXMAX#  BOYHIN#  BOYHAX 

C  **♦ 
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c 

c 


DATA  CHECKS 

CALL  RHNHX(SCALE»  4999.#  lOOOUi*#  ‘SCALc*#  «REAS  MIN*# 
»  *REAS  MAX*) 

CALL  IHNHXCINK#  -1#  2#  'INK*#  ‘NEC  NO*# 

N  *REAS  MAX*) 

RMMHXCXQBS#  XMIN-.OOI#  XHAX+.pOi#  »XQBS*#  *XMIK*# 
RMNMXtYQBSi  YMIN-.ODl#  YHAX+,001#  ’YOBS*#  *YMIN'# 


CALL 

CALL 

CALL 


*XMAX«  ) 
•  YiUX*  ) 


100 


RHNHX(HOBS#  -.CiiCLl# 
•  «REA3  MAX*) 

CALL  IMNMXtlDNTMV#  -I#  2# 
CALL  IHNMX(NPLTHT#  0/  11# 
^  *Pl.THT  DIM*) 

CALL  IHNMXCITGTAV#  -1#  2# 
DO  lot  I  •  1/  NPLTHT 

CALL  RMNHX(PLTHT(n#  - 
K  »N£G  NUMBER'# 

CONTINUE 
CALL  RMNMX(XM1N# 
RNNHX(XMAX# 
RHNMX<YMIN# 
RHNHX(YKAX# 
RHNMX^GRIDR 


10001.#  'HOBS'#  'NEG.  NO.*# 

»IDNTMV»#  «N£G.  NO.*#  *R£AS 
•NO.  OF  PLOT  HTS,*#  » ZERO'# 


MAX*  ) 


•STGTAVS#  *N£G.  NO 


. '  # 


•REAS  MAX') 


.OOOUl#  lOOLl.# 
•REAS  MAX*) 


»NPLTHT{I)*, 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


99999.#  XHAX#  'XHIN*#  'REAS  MIN'#  *XHAX*) 
XMIN#  1000000.#  'XMAX*#  'XHIN*#  *P.£AS  MAX*) 
999999.#YHAX#  'YHIN*#  'REAS  MIN'#  *YHAX*) 
YHIN#ii/uODOua.#  »YHAX*#  *YHIN»#  'REAS  MAX*) 
0.#  1001.#  'GRIDR*#  •ZERO'#  'RcAS  MAX*) 


IHNHX(JVEG#“1>3#*VEG  ADD  IN'# *N£G.  NO*# 'MAX  VAL  »2*) 
RKNMX<TGRID#  C.#  iOOl.#  •TGRID*#  ‘ZERO'#  'REAS  MAX') 
6QZIN2(TGRID# 

IHNMXCNTAPES# 

IHNMXdBQX# 


GRIOR#  *TGRID*#  'GRIDR*) 

0#  13#  *NTAP£S*#  'REAS  MIN'#  'MERGE  LIH') 


IF  (IBOX 
IF  (IBOX 


«£Q. 

.EQ. 


1) 

1) 


CALL 


-1#  2#  'IBOX'#  *NEG.  NO'#  'REAS  MAX') 

CALL  RMNHX(BOXHIN#  XMIN-.OOl#  BOXMAX#  'BOXHIN'# 
•XMIN'#  'BiJXHAX') 

RMNMXtBUXHAX*  BOXHIN#  XHAX+.L^i#  'BOXMAX'# 
•BOXMIN'#  'XKAX') 

RMNMX<BOVMIN#  YHIN-.OOl#  DOYMAX#  'BOYHIN'# 
»YMIN*#  'BDYHAX') 

IF  (IBOX  .EQ.  1)  CALL  RMNMX( DOYMAX#  BOYMIN#  YHAX+.OOl#  'BOYMAX'# 

'BOYHIN'#  'YMAX') 

•ISAVE*#  *NEG.  NO'#  'PEAS  MAX') 


IF  (IBOX  .EQ.  1)  CALL 


2# 


C  4>A4i 

C 

c 

900 

910 

920 

930 

94C 

950 

960 


CALL  iMNHXdSAVE#  -X# 
RETURN 

FORMATS 


FORHATtAlL#  FIO.O#  IlU) 
FORMATdlO#  3F10.C#  110) 
FQRHAKBFIO.O) 
FORMAT(5Flt.C#IiO) 
F0RHAT(iX#A10#F10.0#I10) 
FORMAT (*  ELEVATION  OF  TERRAIN 
FORMATdlO#  110) 

END 


AT  OBSERVER  POSITION  -*#  FlO.l) 


4. 
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SUBRQUTZNe  RMNNX(ZVALUEiZnZK»ZHAX« VAlNAHfHXHNAN^NAXNAN) 
BY  BARBARA  BROOME 

MAKES  SURE  ZHXN  tS  IES$  THAH  ZVALUl  IS  LESS  THAN  ZHAX, 
VALKAH  •  'ZVALUE*#  MINNAH  •  «ZKXN*«  MAXNAH  <»  'ZHAXS 

If  (ZVALUE  tLE*  ZMIN)  WRITE(6«100)  VALKAH/ HINNAH# VALNAH^ 

*  ZVALUE/HINHAH/ZHIN 
IF  (ZVALUE  .GE*  ZHAX)  WRZTE(6/Z00)  VALKAH/ HAXNAH/ VALNAH/ 

«  ZVALUE/ HAXNAH/ ZHAX 

IF  (ZVALUE  .LE.  ZHINl  STOP 
IF  (ZVALUE  .GE*  ZHAX)  STOP 
RETURN 

100  FORHATdX/  A20/ «  SHOULD  BE  GREATER  THAN  •/  AZw/// 

<•  IX/  AZO/*  -  */  FIO.O/  // 

*  IX/  A2D/*  •  •/  flOtO) 

200  FORHATdX/  A20/*  SHOULD  BE  LESS  THAN  •/  A2L/// 

P  IX/  A2p/«  •  •/  F10«>0/  // 

*  IX/  A20/*  •  •/  flO.O) 

END 
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SUBROUTINE  RTYPRC  (ICQDE#  I) 

C  *** 

C  BY  KOHTE  COLEMAN 

C  ***  CHECKS  TO  HAKE  SURE  INPUT  RECORD  IS  THE  R15HT  TYPE 
C  *** 

IF  ((I  iSQ*!  1)  .lANO.  (ICODE  •NE«i  iTAPE  IO«»}  GOTO  100 
IF  ((I  tCQ^L  2)  •ANO«  (ICQOE  •NE»  *FILE  ID*))  GOTO  lOlO 
IF  (ICQDE  .£Q.:  *EQFM  GOTO  200 
IF  (ICODE  .CQ*!  «EQI«)  GOTO  300 

IF  {(I  .GT.t  2)  .AND*  (ICOOE  aNE.*  *DATA> } )  GOTO  100 

IF  (Z  tOT*  2500)  GOTO  400 

RETURM 

C  ***> 

C  ***  THE  FIRST  2  RECORDS  SHOULD  BE  FZLE  ID  AND  TAPE  ZD 
C 

100  URXTE(6>150)  I* ICQDE 

ISO  FORMAT! *1REC0RD  S  I9«  •  SHOULD  NOT  BE  A  *#  A10«  *  RECORD. •) 

STOP 

C 

c  ♦fO  IF  EOF  IS  ENCOUNTERED 
C  >04"^ 

200  WRITE(6#250)  ICUOE^I 

250  FORMAT(//#<  *f  AlO#  »  ENCOUNTERED  AT  RECORD  S  15) 

RETURN 

C  <lli4<4i 

C  IF  EOI  IS  ENCOUNTERED 
C 

300  URITE(6#2501  1C0DE«I 
ST(JP 
C 

C  IF  there  SEEMS  TO  BE  TOO  MANY  RECORDS  ON  THE  TAPE 
C  *** 

400  WRITE! 6# 450 J 

450  rnRMAT(*lARE  YOU  SURE  YOU  SHOULD  HAVE  MORE  THAN  2500  SCANLIMESTM 
STOP 
END 


o  w 


SUBROUriNE  SAVi‘F.:(  ITCTAV) 

>^4<« 

BY  BARBARA  ERQOHE 

SAVL  PRELIHIHARY  FILE  ID  INFO  IN  RECORDS  4001-A013#  IF  ISAVE  •  1. 

COMMON  GRIDR 

OOHKOM  /SCAMS/  MX>  HY 

CQHMOM  .'SEE/  YOSSf  HOBS#  OBTRAN 

COHhOH  /AREA/  XtlRiCN,  YUP.IGM/  XSTOP#  YSTOP 

COMMON  /PLTSTF/  MAPID,  SCALE#  NPLTHTt  PLTHT(10)»  lOBS#  INK 

Ultir.NSinM  MIMOEXtAOlS) 

CALL  OPEMHSUtf  MINDEX#  4015#  0) 

CmLL  FRITHS XOBSi'  1#  4001) 

CALL  WRITNS  (5v  YOBS#  1»  4002) 

CAL^  L'RITHS  i3,  MOBSi  1#  4003) 

CALI  WRUMS  (3*  XORIGNj  1#  4004) 

CALL  WP.ITMS  (3.  YORIGIl#  1»  4005) 

CAlL  L'P.ITMS  Ut  XSTOP#  1#  4006) 

CALI  WR:T0S  (3#  YSTUP#  1#  4007) 

CALL  CRITIIS  (3,  MX#  1#  4006) 

CALL  kif.-THS  <3#  KY#  1#  4009) 

CALL  \;P.j’.TMS  <3#  GI'.IC-R#  1#  4010) 

CAILWr.IVhS  (3i  MAPID#  1#  4011) 

CAI  I  NRITiiS  <3#  IQBS#  1#  4012) 

CALL  X'RITHS  i3#  ITCTAV#  1#  40X3) 

r.M'OFN 

El  ID 
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SUBROUTINE  ST0ATA<Z«J«N««1VE8) 


4i4>4i 

BY  ARTHUR  6R0VES 

STORE  ELEVATIONS  PLUS  VEOETATION  PLUS  INDICATOR  FRACTIONS  IN 
ARRAY  Z  STARTING  AT  Z(4)»i  N  ■  DXNENSION  OF  Z-ARRAY. 

COHHON  /DHACON/  XBUF(700}»ITENP(B0)#ZTEH 
DIHENSIQN  Z(1)#VECCQ0(4) 

INTEGER  SHIFT 

DATA  HSK1/10B/»HSK2/77777B/*VEGCOD/0.0»20.3#10.Z«A.1/ 

I*J 

lUORO-Z 

IGET-l 

CALL  UNPACK(IBUF(IMQRD)«ITEHP»60) 

ZU0R0»ZUaRD4'6 

100  XSIGN»ZTEHP{I6ET).AN0.MSK1 
IVEG-SHIFTCITEHP(IGET)#-^) ^ 

1TEH»ICVT(ITEHPC1GET)«3)  .and*  HSK2 
IF ( ZTEH«EQ«HSK2) RETURN 
IF<ISIGN«NE«0)ITEH— ITEM 
Z(!)«FLaAT(ITEH} 

IF(JVEG«EQ«0)Z(I)-Z(Z)4VE6CODaVEG«>l) 

IFiJVEG.EOa  .AND.l  IVEG«CQ»2}Z(X)-Z(I)^VEGC0D(3) 

I-I+l 

ZF(I«GT.N)RETURN 

IGET-XGET+3 

1F(H0D(IGET«6Q)«NE.1)G0T0  100 
CALL  UNPACK(IBUF(luaRD)fITEHP*60) 

ZUaRO-IWORD4-6 
IGETrl 
GOTO  100 
END 
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FUHCTIon  TVEG(2) 


by  ARTHUR  GROVE3  ^  -tucu  7 

RETURNS  VEGETATION  AHDUNT  ADDED  FOR  THE  GIVEN  Z 


data  QRCHAD  /A,a/»  urban  /10.2/7  FOREST  /20#3/ 
DPZ  •  Z  -  AINT(Z) 


TVEG  ■  0 

IF  (HPZ  .GT.  .05)  TVCG  -  QRCHAD 

IF  (DPZ  .GT.  .15)  TVE6  -  URBAN 

IF  (DPZ  .GT.  .25)  TVE6  -  FOREST 

RETURN 
END 
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SUBROUTINE  XPOSXT(X«Y*ICAREI 

#4«i> 

BY  ARTHUR  BROVES 

THE  PURPOSE  OP  THIS  SUBROUTINE  IS  TO  ATTEMPT  TO  HOVE  AN  OBSERVER 
AT  (X*Y)  NHXCH  is  XN  VEOETATION/URBAN  TO  A  NEARBY  POSITION  WHICH 
IS  IN  THE  OPEN*!  TO  HAVE  ST  DO  THISv  CALL  THE  SUBROUTINE  WITH 
ICARE-0«(  XT  then  LOOKS  AT  THE  ELEVATIONS  AT  THE  FOUR  CORNERS  OF 
THE  SQUARE  OF  ELEVATIONS  (X*Y)  IS  IN«  IF  NONE  OF  THESE  ARE 
WOQDED/URBAN#  IT  RETURNS  THE  ORIGINAL  (X#Y)«.  IF  ALL  FOUR  CORNERS 
ARE  WQOOEO/URDAN*  OR  IF  EXACTLY  ONE  PAIR  OF  OPPOSITE  CORNERS  ARE 
UOOOEO/URBAN*  XT  PRINTS  AN  APPROPRIATE  STATEMENT  AND  STOPS  THE 
RUN#  MAKING  NO  ATTEMPT  TO  FIND  AN  OPEN  NEARBY  LOCATION.  IF  ONE# 
TWO  (ADJACENT)  OR  THREE  CORNERS  ARE  HOODED/URBAN#  IT  MOVES  THE 
OBSERVER  THE  LEAST  POSSIBLE  DISTANCE  TO  AN  ADJACENT  SQUARE  NOT 
KNOWN  TO  HAVE  ANY  WOODEO/UXBAN  CORNERS**  IF  THIS  NEW  SQUARE  XN 
FACT  HAS  NO  WOODED/URBAH  CORNERS*  XT  RETURNS  AS  (X#Y)  THE  NEW 
ADJUSTED  OBSERVER  LOCATION  AND  PRINTS  THE  NEW  LOCATION*  IF  THIS 
NEW  SQUARE  DOES  HAVE  ONE  OR  HORE  WOODED/URBAN  CORNERS#  IT  PRINTS 
THAT  IT  HADE  AN  UNSUCCESSFUL  ATTENPr  TO  NOVE  THE  OBSERVER  AND 
THEN  STOPS  THE  RUN*  IF  THE  SUBROUTINE  IS  CALLED  WITH  XCARE-1# 
THIS  ENTIRE  PROCEDURE  IS  BYPASSED  AND  (X#Y)  IS  RETURNED  UNCHANGED 
WHETHER  OR  NOT  VEGETATION  OR  URBAN  ARE  PRESENT. 

*4141 

DIMENSION  T(A#A) 

COMMON  GRXDR 
COMMON  /SCANS/  NX#NY 

COMMON  /AREA/  XQRX6H#  YQRI6N#XSTOP« YSTOP 
COMMON  /BIG/  SCAN(4000)  ^ 

LOGICAL  Ll#L2fL3#L4 


1  FORMAT(«  ELEVATIONS  REQUIRED  FOR  XPOZXT  NOT  PROVIDED  BY  MERGIT') 

Z  FOKMATC  OBSERVER  MOVED  TO  (•#F10*a#  *#*#F10.1#«)  •) 

3  FCRMAT(«  OBSERVER  ENTIRELY  SURROUNDED  BY  VEGETATION  OR  URBAN  AND  C 
*  OUI.D  NOT  BE  MOVED*) 

4  FORMATC  EXACTLY  ONE  PAIR  OF  OPPOSITE  CORNERS  IN  VEGETATION  OR  URB 
«  AN  -  OBSERVER  COULD  HOT  BE  MOVED') 

5  FORMATt*  HADE  ONE  UNSUCCESSFUL  ATTEMPT  TO  MOVE  OBSERVER*) 

*4i4‘ 

If  NO  ArrtHPT  TO  MOVE  THE  OBSERVER  IS  TO  BE  MADE#  RETURN. 

*4'* 

XFdCARE.EQ.DRETURM 

4<4'4« 

IDENTIFY  THE  LOCATION  OF  THE  4X4  SQUARE  OF  ELEVATIONS  SURROUNDING 
THE  OBSERVER. 

**41 

ILoLOCATECX* XORIGH^GRIDR) 

Jl»LOCATE(Y#VORISN#GRXDR) 

XORIG  «  XORIGrS  ❖  CRIDR  *  FLOAT? XL-2) 

YORIG  ■  YORICN  ♦  GR2DR  ♦  FLOAT (JL'*2J 

IF  ALL  16  OF  THESE  ELEVATIONS  HAVE  NOT  BEEN  PROVIDED  BY  THE 
MERCIT  SUBROUTINE#  PRINT  THAT  FACT  AND  STOP  THE  RUM. 


C 


IF(XL»l.GEoa.AN0*lL*-2.LE*NX*AND»UL-l«G£.l*AND.JL«2.L£«NY)GaT0  100 
WRZT£(6#1} 

STOP 


oooo  nc>r>r>rt  r>ooooo  r>oor»oo  r»ooo 


C  FILL  THE  4X4  ARRAY  WITH  ELEVATIONS, 

C  4<4i* 

100  DO  120  z«l,4 

CALL  REA0nS(2#$CAN<l)#HY«£L«'Z'^2) 

DO  110  J»l«4 

Ta#j)«S&AN(JL4-J-2) 

110  CONTINUE 

120  COHTINUE 

C  >*** 

C  CLEAR  THE  •SCAN*  ARRAY,  JUST  IN  CASE* 
C  *** 

OD  130  Z«l«4000 
SCAN(Z)»0* 

130  CONTINUE 


SET  •NQVE0<  TO  ZERO  TO  INDICATE  THAT  (X^Y)  IS  THE  INPUT 
OBSERVER  POSITION  AND  NOT  AN  ADJUSTED  POSITION* 

««« 

HQVED-0 

#4»li 

OETERHINE  WHICH  CORNERS  OF  THE  2X2  SQUARE  CONTAINING  THE  CURRENT 
OBSERVER  POSITION  ARE  IN  VEGETATION  OR  URBAN*  THE  APPROPRIATE 
LOGICAL  VARIABLE  < LI,  L2,  L3  OR  L4)  WILL  BE  HADE  •TRUE*  IF  THE 
CORRESPONDING  CORNER  IS  IN  VEGETATION  OR  URBAN. 

140  I-L0CATE(X,X0RI6,GR10R) 

J«LOCATE(Y,YQRIG,GRIDR) 

XSCAN  -  XORIG  ♦  GRIDR  ♦  FL0AT<I-1) 

YSCAN  -  YORIG  ♦  GRIDR  ♦  FLQAT<J-1) 

LI  -  AIMT(10*PFRACT(T(I,J)*-*00001)}  .NE*  0* 

L2  ••  AXNT(10*P*>FRACT(T(Itl,J)nD0001))  .NE*  0* 

L3*  AINT(10*t»FRACT(T(I>l,JnU*D0001)}  .NE*  0* 

L4  -  AlNT(10*«FRACt(Ta,J4l)nOOOOin  .NE.  0* 

♦♦♦ 

IF  NO  CORNERS  ARE  IN  VEGETATION  OR  URBAN#  THE  OBSERVER  IS  IN  THE 
OPEN*  IF  THIS  IS  THE  ORIGINAL  OBSERVER  POSITION#  RETURN  TO  THE 
CALLING  PROGRAM.  IF  THIS  IS  AN  ADJUSTED  POSITION#  WRITE  THAT  FACT 
AND  THEN  RETURN 

IF  (Ll.CR*L2.nR.L3.lQR«L4)  GOTO  150 

1F(K0VED*EQ.1)WRITE(6,2)X«Y 

RETURN 

**41 


IF  HQVED-1  INDICATING  THAT  THE  OBSERVER  HAS  ALREADY  BEEN  MOVED 
TO  A  NEW  LOCATION  AND  IS  STILL  IN  VEGETATION  DR  URBAN#  WRITE 
THIS  FACT  AND  STOP  THE  RUN.' 

*** 

150  IF(HQVED.EQ.0)G0T0  160 
WRITE(6,5) 

STOP 

**♦ 

IF  OBSERVER  IS  ENTIRELY  SURROUNDED  BY  VEGETATION  OR  URBAN# 
WRITE  THIS  FACT  AND  STOP  THE  RUN. 

*** 


160  IF(*N0T.'(Ll.AND.12*AND.L3.AND.t4nG0T0  170 
WRITE(6,3) 

STOP 
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c 

C  IF  ONE  PAIR  OF  OPPOSITE  CORNERS  ARE  IN  VEGETATION  OR  URBAN#  WITH 

C  THE  OTHER  PAIR  IN  THE  OPEN#  WRITE  THIS  FACT  AND  STOP  THE  RUN, 

C 

170  IF(,NOT*iai,ANO«L3*AND.iNOT«L2*ANO««NOT«U  .OR.  L2«AND«I.A*AND* 

*  •NQT«L1«AND,WN0T«L3))G0T0  ISO 
WRITE(6#4) 

STOP 


C  *** 

C  THE  OBSERVER  IS  GOING  TO  BE  MOVED  TO  A  NEW  LOCATION*  AT  THIS 
C  POINT  HAKE  'MOVED'  EQUAL  1  SO  THAT  WHEN  CONTROL  RETURNS  TO 

C  STATMENT  140  TO  EXAMINE  THE  NEW  LOCATION#  IT  WILL  RECOGNIZE  IT 

C  AS  A  MEU  LOCATION.i 

C  *** 

180  HOVED-1 

C  ♦♦♦ 

C  HOVE  THE  OBSERVER  IF  CORNER  NUMBER  1  IS  THE  ONLY  CORNER 
C  IN  VEGETATION  OR  URBAN, 

C  ♦♦♦ 

IF(,MOT..{  Ll,AHD,*NaT,L2,AH0*.N0T*L3.ANP,.NQT*L4))GnT0  190 

Dl-XSCAN  ♦  GRIDR  -  X 

nZ-YSCAN  ♦  GRIDR  -  Y 

IF(D1,LT,02)X-XSCAN  *  GRIPR  ♦  ,01 

IF(Dl,CE,D2)Y-YSCAN  *  GRIDR  +  ,01 

GOTO  140 

C  *♦4 

C  HOVE  THE  OBSERVER  IF  CORNER  NUMBER  2  IS  THE  ONLY  CORNER 
C  III  VEGETATION  OR  URBAN* 

C  *** 

190  ir(,NQT,l(*NUT*Ll*AND*L2*AM0**N0T*L3»AND**N0T*L4)  )60T0  200 
Dl-X  -  XSCAM 
nZ-YSCAN  ♦  GRIDR  -  Y 
IF(D1.LT*02)X»XSCAN  -  *01 
IF(D1.CE*02)Y-Y5CAN  +  GRIDR  ♦  ,01 
GOTO  140 

C  44* 

C  MOVE  THE  OBSERVER  IF  CORNER  NUMBER  3  IS  THE  ONLY  CORNER 
C  III  VEGETATION  OR  URBAN, 

C  444 

200  ir(,NOT,l(,NOT,Ll*AND,,HOT,L2,AND,L3,AND,,KaT,L4)  )GOTO  210 
Ul-X  -  XSCAN 
02 -Y  -  YSCAN 

IF <D1,LT,02JX«XSCAN  -  ,01 
IF(D1.CE,02)Y-YSCAN  -  ,01 
GOTO  140 


C  444 

C  MOVE  THE  OBSERVER  IF  CORNER  NUMBER  4  IS  THE  ONLY  CORNkR 
C  IH  VEGETATION  OR  URBAN* 

C  444 

210  IF{,NaT,l{*NDT,Ll*AN0,,NnT,L2,AND*jNGT,L3,AM0,l4)  )GOTO  220 
01-XSCAN  ♦  GRIDR  -  X 
D2-Y  -  YSCAN 

ir(Pl,LT,D2)X-XGCAN  ♦  GP.IPR  +  ,01 
IF(D1,CE«D2)Y-YSCAH  -  ,01 
GOTO  140 

C  444 

C  MOVE  THE  OBSERVER  IF  CORNERS  NUMBER  1  E  2  ARE  THt  ONLY 
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C  CORNERS  ZN  VEGETATION  OR  URBAN«i 

C 

220  ir<.NOT*i(Ll«AND«i2>ANO«.NOT«U«AND«.NOT«L4))GOTO  230 
Y-YSCAN  *  GRXDR  *  *01 
GOTO  140 

♦♦♦ 

HOVE  THE  OBSERVER  IF  CORNERS  NUHBER  1  C  4  ARE  THE  ONLY 
CORNERS  IN  VEGETATION  OR  URBAN.i 

230  IFUNQT»l(Ll«AND«l»NOT,L2uAHD*4NOTeL3»AND»L4))60TO  240 
X-XSCAN  ♦  GRIDR  *  .01 
GOTO  140 

*** 

HOVE  THE  OBSERVER  IF  CORNERS  NUHBER  2  C  3  ARE  THE  ONLY 
CORNERS  IN  VEGETATION  OR  URBAH.i 

4i4t4i 

240  IF(.N0T«IC«NQT«Ll«AND«L2tAN0«t.3«AND»«N0T*L4))G0T0  250 
X-XSCAN  -  *01 
GOTO  140 

HOVE  THE  OBSERVER  IF  CORNERS  NUMBER  3  C  4  ARE  THE  ONLY 
CORNERS  IN  VEGETATION  OR  URBAN, I 

250  IF(«N0T«l(aHQT«il«AND««N0T,t2«AHD»L3*'AND«L4))G0Tn  260 
Y-YSCAN  -  ,101 
GOTO  140 

HOVE  THE  OBSERVER  IF  CORNER  4  IS  THE  ONLY  ONE  IN  THE  OPEN, 

260  IF(,NOT,l(Ll,AND,L2,AHO,L3.ANO,i,NaT,L4)IGOTn  270 
X-XSCAH  -  ,01 
Y-YSCAN  ♦  GRIDR  ♦  ,01 
GOTO  140 

MOVE  THE  OBSERVER  IF  CORNER  3  IS  THE  ONLY  ONE  IN  THE  OPEN, 

4<4>4> 

270  IF(,NQT«l(Ll,AND,C2«AHl>,»N0T,<.3,AHDoL4nG0Ta  280 
X-XSCAM  ♦  GRIDR  ♦  ,01 
Y-YSCAN  GRIDR  +  ,01 
GOTO  140 

4<4i4i 

HOVE  THE  OBSERVER  IF  CORNER  2  IS  THE  ONLY  ONE  IN  THE  OPEN, 

*** 

280  IF(,NaT,  (Ll,ANn,i,NGT,L2,AND,i3,AND,L4))G0Tn  290 
X-XSCAN  ♦  GRIDR  ,01 
Y-YSCAN  -  ,01 
GOTO  140 

♦♦♦ 

MOVE  THE  OBSERVER  IF  CORNER  1  IS  THE  ONLY  ONE  IN  THE  OPEN, 
THIS  IS  THE  ONLY  CASE  LEFT, 

290  X-XSCAM  -  ,01 
Y-YSCAN  “  ,01 
GOTO  140 
END 
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rUHCTlUi^  ’rRiriE(XTGT#  YTCT#  XTGT>  ZOLD#  ZLAT2R) 

BY  BARBARA  BROCIHE 

IF  ZPRIHE  IS  GREATER  THAN  THE  TARGET  ELEVATION*  PREVIOUS  TERRAIN 
HUST  BE  BLOCKING  THE  VIEW*)  TGT  VISIBLE  IFF  HOUHI  nlE*  0. 

ZOLD  array  COMTAIMIHG  ELEVATIONS  FOR  PREVIOUS  SCAN  (TRUMPED  UP), 
ZLATER  WILL  BE  ZOLD  FOR  MEXT  SCAN  LINE 

X,Y  COORDINATES  OF  INTERSECT  OF  0BS-T6T  LINE  WITH  ZOLD  SCAN. 

Z  interpolated  elevation  corresponding  TO  POINT  (X*Y) 

RE  RADIUS  OF  THE  EARTH  IN  HETERS 
ZPRIHE  PROJECTION  OF  Z  TO  TGT  POINT 

DIMENSION  ZULDiDi  ZLATER'*) 

COMMON  CRXPR 

COMMON  /SEE/  XOBS*  YOBS*  HOBS#  OBTRAN 

LOGICAL  QUAOlt'#  QUADlD#  QUAD2U#  QUAD20#  QUA03U#  QUA03D* 

*  QUA04U#  QUA04D 

RE  -  6376323, 

DX  -  KTGT  -  XQBS 
OY  «  YTGT  -  YOBS 

♦•ii# 

DEFINE  X#  y»  Z  AFTER  DETERMINING  WHICH  PORTION  OF  WHICH  QUADRANT 
YOUR  TARGET  IS  IN  (3  POS IBILITIES) 

QUADIU  •  (OX  «GE;.  0«)  .AND*  (OY  tGE.i  0, )  •AND.  (DX  .LT*  DY) 

IF (.NOT,.  QUAD) U) GO  TO  150 
X  -  OX*(DY-ORIOk)/OY  ♦  XDDS 

Y  •  YTGT  -  GRIDR 

Z  •  (l-0X/0n*ZLATER(IT0T-l)  ♦  ( DX/OY)»ZOLDi  ITGT-1) 

GO  TO  500 

C  **♦ 

ISO  QUADlD  -  (OX  «GE«  0.)  .AND.)  (DY  .GE,.  0>)  ,AND,  (OX  ,GE,  OY) 
IF(.Nar,i  QUADiDJGO  TO  200 
X  •  XTGT  -  GRIOR 

Y  -  nYPC^X-OklOPJ/OX  YOBS 

Z  •  (1-C'\'OX)*ZOLO(ITGTJ  ♦  ( DY/DX) *20LD ( ITGT-1 ) 

GO  TO  50 <• 

C  ♦** 

200  QUAnaU  »  (OX  .LT.  O.)  .ANO.i  (DY  .GE,  0.)  .AND,  C-OX  .LTq  DY) 
IF(oNDT,'  CUAD2ii)G0  TO  250 
X  -  0X*(DY-CRIDRJ/DY  ♦  XOBS 

Y  «  YTGT  -  GRIDR 

Z  -  (l+DX/DY)*ZLATEP.(i:TGT-l)  ’•  ( DX/DY)*ZOLD(  ITGT-l) 

GO  TO  500 
C 

250  QUADZl'  -  (DX  .LT*  C. )  .AND,.  (DV  ,G.E,  C, )  .AND,  (-DX  «GE.  DYI 
IF(*NOT,i  QUAD2D)Ga  TO  300 
X  -  XTGT  ♦  GRIDR 

Y  ■  DY*(OX+GRIDRVDX  ♦  YOBS 

Z  •  (l<-DY/DX:*Zr)LD(ITGT)  -  (DY/OX)*ZOED(  ITGT-1) 

GO  TO  500 
C 

300  Q»JAD3U  •  50X  .LT.  0,)  •  AND,  (DY  ,LT?  0.)  .AND,  (OX  ,LT.  DY) 
IF(.N0T,1  QUA03U)G0  TO  350 
X  •  XTGT  +  GRIOR 

Y  ■  OY*(DX+GRIDR)/DX  +  YOBS 
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2  •  (1-0Y/DX)*Z0LD(ITGT)  ♦  tOY/OX)*ZOLO(nGT+ll 
GO  TO  500 

C  ♦♦♦ 

350  QUAD3D  •  (OX  •(.!•  0«J  •AKO.i  (DY  •LTti  0«)  .AND*  (DX  *GE»  OY) 
IF(«NOT.l  QUAD3D)GQ  TO  400 
X  -  DX4<(DY*GRZDR)/DY  4-  XQBS 

Y  ■  YTGT  ♦  6R1DR 

Z  -  (l-DX/DY)*ZLATEP.{ITGT^l)  ♦  (  OX/DY)  >«<ZaiD(  ITGT  +  1) 

GO  TO  500 

C  ♦** 

400  QUAD4U  »  (DX  .CE*  0.)  .AMD*.  (OV  «LT«  0*}  «AKD»  (DX  *GT«  -DY) 
IFi.NOT.i  QUAD4U)GQ  TO  450 
X  -  XTGT  -  GRIDP. 

Y  ■  DY*(DX-GRIDR)/DX  +  YOBS 

Z  ■  (1«-DY/0X)*Z0LD(ITGT)  -  (DY/DX) ♦ZDLD(  ITGTn ) 

GO  TO  500 
C  ***  QUA040 

450  X  ■  DX*(DY+GRIDR)/DY  ♦  XOBS 

Y  -  YTGT  ♦  GRIOR 

Z  «  {l*DX/DY)4ZLATER(ITGT+l)  -  ( DX/DY) ♦ZOLD{ ITGT+l) 

C  ♦** 

500  CONTIMUE 

C  *** 

C  CALCULATE  THE  PRQJECTIUH  OF  2  ONTO  THE  TARGET  POSITION 
C 

RZ  •  SQRT((X-X0BS)**2  +  (Y-Y0BSi**2) 

RT  -  SQRT<(XT6T-X0BS)*42  ♦  (YTGT-YOBS) **2) 

B  -  RE  +  HOBS  *  OBTRAN 

A  «  ({RE+Z)'>caS<RZ/RE)-RE-HOBS“OBTRAH)  /  ( (RE+2)  ♦SItKRZ/RE ) ) 
XPRIME  ■  B  /  (COT<RT/RE)«A) 

YPRIME  •  COT{PT/RE)*B  /  (CQT(RT/RE)-A) 

ZPRIIJE  •  SQRT(XPRIME**2  ♦  YPRIHE**2>  -  RE 

RETURN 

END 
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requires  1-  NPLTHT  -  10. 

ITGTAV: Indicates  whether  target  height 

is  the  target's  height  above  terrain 
only  (IT6TAV=0)  or  above  terrain 
and  any  vegetation/urban  (ITGTAV=1). 


CARD  :  INPUT  i  FORMAT  ‘=  DESCRIPTION 


RDNCHK 


READ  AND 
\CHECK  INPUT, 


FRSTiM 


FOR  2 
SCANS  NEAR-' 
/EST  OBSERVER 
CALCLMTE 
HOWHI's  & 
HORIZON'S 


MERGIT 

PUT  BfV- 
ATIONS  IN 
DIRECT  ACC-j 
ESS  FILE 


APPENDIX  A3 
SEEFAR  FLOWCHART 
XPOSIT  aVOBS 


IP  neb:)ED, 
AAOVE  OBSER¬ 
VER  OUT  OF 
.OBSTRUCTION 


CALQXATE 

06SERVB)'S 

ELEVATION 


PLTPRE 

PUT  PRE 
LIMINARY 
INFO  IN 
PLOT  FILE 


SAVPRE 

IF  NEHDe), 
SAVE  PRELIM¬ 
INARY  INFO  IN/ 
HOLD  FILE 


PLTLN 


PUT  LOS 
-H  RESULTS  IN] 
iPlDT  FILE 

/if  needed,) 

/UPDATE  NU- 
/MERATOR  & 
iDENOMINATC 
\FOR  VISIBLE 
FRACTION 


WRITMS 

IF  NEEDED,  \ 
-H  PUT  HOWHI's 
,  IN  HOLD  FILE  / 
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Appendix  B1 

Calculation  of  Critical  Horizon  Coordinates 

In  determining  whether  a  target  position  is  in- view  or  out-of- 
view  the  first  step  is  to  consider  the  effect  of  intervening  terrain. 

This  is  done  by  finding  where  the  observer- target  line  intersects  the 
closest  line  of  known  horizon  values  (as  illustrated  in  Figure  Bl-1). 

This  intersection  location  is  referred  to  as  (X,Y).  To  find  the  horizon 
value,  Z,  at  (X,Y),  one  need  only  interpolate  between  the  closest  hori¬ 
zon  values  on  the  horizon  scan,  pictured  as  Z^,  I2' 

The  problem  of  solving  for  X,Y,  and  z  breaks  into  eight  cases 
(illustrated  in  Figure  Bl-2).  These  cases  are  determined  by  the  relation¬ 
ship  between  the  coordinates  of  the  target  and  those  of  the  observer. 

Considering  a  pair  of  axes  with  origin  at  the  observer  position, 
each  quadrant  has  been  divided  into  two  sectors.  Targets  in  one  sector 
have  an  observer- target  line  slope  whose  absolute  value  is  less  than  or 
equal  to  one.  In  the  other  sector  the  absolute  value  of  this  slope  is 
greater  than  one.  The  reason  for  this  division  will  be  explained  by 
the  following  two  examples. 

Selecting  quadrant  1,  X,  Y,  and  Z  will  be  determined  for  the 
area  where  XTGT  ^  XOBS,  YTGT  ^  YOBS  and  (XTGT  -  XOBS)  >  (YTGT  -  YOBS). 

(See  Figure  Bl-3. ) 


Letting  GRID  be  the  grid  of  the  map  data  X  =  XTGT  -  GRID,  since 
the  horizon  is  one  scan  away.  Knowing  X,  we  can  solve  for  Y  from  the 
equation  of  the  observer-target  line: 

Y  =  (YTGT- YOBS)  (XTGT-XOBS-GRID)  +  YOBS. 

(XTGT- XOBS) 

It  remains  to  solve  for  Z.  Consider  the  target  to  be  the  ITGT^^ 
position  on  a  scan  and  ZOLD  to  be  an  array  containing  old  horizon  values. 


Let 


DX  =  XTGT-XOBS  and  DY  =  YTGT-YOBS. 


Then  linearly  interpolating, 
Z  =  (1  - 


^)  ZOLD  (ITGT)  +  4^)  ZOLD  (ITGT-1) 


Now  consider  a  target  in  the  upper  portion  of  quadrant  1,  where 
XTGT  -  XOBS,  TYGT  ^  YOBS  and  (XTGT-XOBS)  <  (YTGT  -  YOBS).  See  Figure  Bl-4.) 
In  this  example  the  closest  line  of  known  horizon  values  is  formed  by  ZOLD 
(ITGT-1)  and  ZLATER  (ITGT-1)  rather  than  by  two  old  horizon  values.  In  this 
case 


Y  =  YTGT  -  GRID 
X  =  DX  (DY-GRID)  +  XOBS 


(1 


DY 


-  DX)  ZLATER  (ITGT-1)  +  -(w)  ^OLD  (lTGT-1) 
DY 


Where  DX,  DY  and  GRID  are  as  defined  in  the  previous  example.  Similarly, 
equations  for  X,  Y,  and  Z  may  be  derived  for  each  of  the  eight  sectors, 
as  indicated  in  Figure  Bl-5. 
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HORIZON  SCAN 


Figure  Bl-1 

Locating  Intersection  of  Observer  -  target  line  with 
Closest  Horizon  Values 


XTGT  <  XOBS 


XTGT  >  XOBS 


YTGT  >  YOBS 


YTGT  >  YOBS 


(XTGT-XOBS)<  (YTGT-YOBS)  (XTGT-XOBS)  <1  CYTGT-YOBS) 


XTGT  XOBS 

YTGT  ^  YOBS  ^  ^ 

-(XTGT-XOBS)  >  (YTGT-YOBS)^ 


XTGT  <  XOBS 
Yl’GT  <  YOBS  ^ 

(XTGT-XOBS) <  / 

CYTGT-YOBS)  / 


XTGT^  XOBS 
YTGT  ^  YOBS 

(XTGT-XOBS)  >  (YTGT-YOBS) 


XTGT  <  XOBS 
YTGT  <  YOBS 


s  CXOBS, 

YOBS) 

XTGT  >  XOBS 

s 

YTGT  <  YOBS 

\ 

\ 

(XTGT-XOBS)^-  (YTGT-YOBS 

\ 

S 

\ 

XTGT 

^  XOBS 

\ 

YTGT 

<  YOBS 

\ 

(XTGT-XOBS)  <  - 

\ 

(YTGT-YOBS)  v 

Figure  Bl-2 

Eight  Cases  for  Computing  X,Y,Z 


ZOLD  SCAN 
(CURRENT  HORIZON) 


Figure  Bl-3 

Solving  for  X,Y,Z  when  XTGT  ^  XOBS,  YTGT  >  YOBS  AND 
(XTGT  -  XOBS)  >  (YTGT  -  YOBS) 


[illlil 


SCAinI 

lORIZO] 


(XOBS.YOBS) 


Flgur*  SI'S 

for  X.Y  'j  for  uho  ITGX^^  Targot  Position  on  s  Scan  (Eight  Possibla  Casas^ 
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2C 

OX<0 

DYiO 

-DXiDY 

XTGT  +  GRID 

^1*5^0LD(ITGT) 

-^ZOLDCITCT-l) 

3U 

DX<0 

DY<0 

DX<DY 

- 1 

XTGT  +  GRID 

P£iPX^H£).Y0M 

(l-J^ZOLD  CITGT) 

♦  ^Z0ID(ITGT+1) 

■ 

3D 

DX<0 

DY<0 

DXiDY 

DXCDY+GRID) 

“TT - *vais 

YTGT  +  GRID 

^-^^UTBRCITGT+1) 

+  ^^0LDCITGT+1) 

4U 

DXiO 

D”<0 

DXI-DY 

XTGT-GRID 

DY(DX-GRID.),yp^ 

(a*pyyULP(ITGT) 

-^ZOLDCITGT+1) 

4D 

DXiO 

DY«0 

DX<-DY 

L_ 

2xaYi2W22.xoL.s 

YTGT  +  GRID 

^1+^2UTER(ITG1+1) 

-^^ZOLD(ITCT+l) 
j _  _ _ J 

(XTGT,  YTCT)  ■  thB  TARGET  COORDINATES 
rxOBS,  YOBS)  I  'IHB  OSSisRVCR  COORDINATES 
GRID  -  MAP  GRID 
DX  »  UGT  -  XOBS 
DY  *  YTGT  -  YOBS 

Z,  £)  m  array  holding  horizon  RBPRgDBNTING  PREVIOUS  ELEVATION 
ILATER  -  Flfi'jRB  HORIZON  VALUES  CALCULATED  POR  CURREKT  SCAN 
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Appendix  82 

In  Appendix  81,  a  method  was  shown  for  calculation  of  x,  y,  and  z,  where 
(x.  y)  is  the  intersection  of  the  observer-target  line  with  the  last  running 
horizon  line  and  z  is  the  horizon  elevation  at  the  position  (x,  y).  This 
Appendix  describes  the  method  for  projecting  the  horizon  elevation  to  the 
target  position.  It  is  at  this  time  that  earth  curvature  is  introduced  into 
the  line  of  sight  calculations. 

Letting  *  the  radius  of  the  earth 

1=  the  elevation  of  the  observer 
(Jd 

H  =*  the  height  of  the  observer 


R  =  the  range  to  the  horizon  line  from  the  observer 
Z 

R  =  the  range  to  the  target  from  the  observer 
the  situation  might  be  drawn  as  follows: 


(0,Rg  +  +  EQg 


(Rg+Z,  II/2  -R^/Rgl 

y 


(V.  w) 


I 


Now 

Z  = 

so  if  we  can  determine  values  V  and  W,  we  can  find  Z  . 

But  (V,W)  is  simply  the  intersection  of  lines  (0,0)  -  (V,W) 
and  (0.  +  E^g)  -  (V.  W). 

The  equation  for  (0,0)  -  (V,W)  is 
Y'  =  AX'  +  B 

where 

B  =  0 

A  =  tan  (n/2-  (''t/R^))  =  cot  (’^T/R^) 

And  the  equation  for  line  (0,Rg  +  HQg+  E^g)  -(V,  W)  is 
Y'  =  A'X'+B' 

where 

"^£'*’*^06''’  ^OB 

and  (R^  Z)  cos  (^^/R  ^  ~  ^OB  ^OB^ 

A'  =  - i - 

(Rg  +  Z)  sin  (R^/Rg) 

Solving  these  two  line  equations  simultaneously,  then,  for  V  and  W 

B' 

^  "  cot  ('^T/Rg)  -  A' 
cot  ('^T/Rg)  B' 

Vj  - ■  ■o  ^ - 

cot  (VRg)  -  A' 
with  A'  and  B'  as  defined  above. 
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APPENDIX  C  1 

The  Time  Complexity  for  an  NXN  Map  Using  the  LOSMAP  Algorithm 


Consider  a  Map  with  N  scan  lines,  and  N  points  per  scan  line  and  an  ob¬ 
server  at  a  point  such  than  n  scan  lines  are  to  his  "left"  and  m  points 
on  each  scan  are  below  him. 

In  example  above,  N  =  7 

n  =  2 

m  =  3 


If  a  calculation  has  to  be  made  each  time  a  profile  crosses  either  a 
vertical  or  horizontal  scan  line,  how  many  calculations  have  to  be  made 
to  compute  a  LOSMAP? 


There  are  N-n  scan  lines  to  the  right  of  the  observer.  The  first 
one  (the  one  immediately  to  the  observer's  right)  requires  N  crossings  of 
vertical  scan  lines  (one  for  each  point  on  the  scan);  the  second  requires 
2N  crossings  of  vertical  scan  lines,  etc.  for  a  total  number  of  crossings 
of 

N  +  2N  +  3N  +  ....  +  (N-n)  N 
=  N  0  +  2  +  3+  ....+  (N-n)} 


N  (N-n;i(N-n+l) 


There  are  n  scan  lines  to  the  left  of  the  observer. 

The  total  number  of  crossings  of  these  is 

N  +  2N  +  ....+  n  N 
=  N  (1+2+. . . .+n) 

=  N^n(n+1)  j 
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I  i..i  . . .  .  i  -j  ' 

,, _ 

■5 

I'l 

j 

I 

Therefore  the  total  number  of  crossings  of  vertical  scan  lines  is 

"  Y  ((N-n)  (N-n+1)  +  n(n+l)} 

n  =  ^  {N^-nN+N-nN+n^-n+n^+n> 

^  2 

=  N  {N2+N+2n^-2nN} 

2 

•  =  N  {N^^-N+2n  (n-N)} 

I 

Now  2n  (n-N)  is  most  negative  when  n  =  N/o 

^  > 

Since  f  (n)  =  2n^  -  2nN 
f (n)  =  4n  -  2N 

Setting  4n  -  2N  =  0,  n  =^/2 

So  f(n)  has  its  minimum  or  maximum  at  N/2 

But  f"(N/2)  =  4>0  so  f(n)  is  minimal  when  n  =  N/2. 

So  n^.  >  ^V2  {N^  +  N  +  2(^/2)  (-^/2)} 

>  ^/2  {N^/2  +  N} 

>  N^/4  ^  N^/2.  order 

We  have  only  looked  at  vertical  scan  lines.  There  are  as  many  again 
crossings  of  horizontal  scan  lines,  (similar  analysis  using  m  in  place 
of  n).  This  essential ly^doubles  the  number  of  calculations,  but  it 
still  remains  of  order  N"'  for  a  map  with  minimum  crossings;  that  is,  of 
order  >M^. 

There  are  always  less  than  2N  crossings  (horizontal  and  vertical)  for 
each  of  the  N^  target  positons,  so  the  order  is  <N^. 

3 

Therefore  the  order  is  N  for  the  LOSMAF  algorithm. 
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